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The operation of a fosstl-fueled thermal 
system (20) is quantified by obtaining a ref- 
erence Kiel chemisoy before on-line opeiatton 
(270). and thereafter operating on-line. In 
on-line operation (270). a set of measurable 
operating parameters is measured, including at 
least effluent concentrations of oxygen and car- 
bon dioxide, and opdonally the concentration 
of effluent water and the concentration of ef- 
fluent sulfur oxide. An indicated Aii/Fuel ratio 
is obtained, as are the ambient concentration 
of oxygen, and air pieheater (36) leakage (29) 
and dilution factors. Hie fuel ash (Eq. 29) 
and fuel water are calculated, and the complete 
As-Fired fiiel chemistry is calculated. From 
the complete As-Fired fuel (Eq. (13)) chem- 
istry, the pertinent systems parameters such as 
reference fuel heating value, boiler efficiency 
(32XEq. 4(20). system efficiency, fiiel flow 
race (Eq. 4(21)). U)tal effiuem flow rate (20). 
individual effluent flow rates (292), and indi- 
vidual emission rates (292) are determined in a 
fully consistent manner. 
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OPERATION OF A FOSSIL-FUEUED THERMAL SYSTEM 



5 This invention relates to a fossil-fired power or steam generation thermal 

system, and more particularly, to a method for determining its fuel chemistry, Tiiel 
heating value, fuel flow, and thermal performance from its basic operating 
parameters. 

The importance of accurately determining thermal efficiency is critical to 

10 any thennal system. If practical day-to-day improvements in efficiency are to be 
made, and/or problems in thermally degraded ^equipment are to be found and 
corrected, then accuracy in determining thermal-efficiency is a necessity. The 
tracking of the efficiency of any thermal system lies fundamentally in measuring 
the useful output and the total in-flow of fiiel. The mefiil output from a fossil 

IS fueled system includes the generation of electrical and/or mechanical power 
and/or the production of a heated working fluid such as^team. 

The measuring of the useful output of themiai systems is highly developed 
and involves the direct measurement of electrical output and/or mechanical drives 
and/or thermal energy flow. Measuring electrical and mechanical power is well 

20 established. Measuring thermal energy flow involves direct measurement of the 
inlet and outlet pressure, temperature and/or quality of a fluid being heated by the 
combustion gases, as well as measurement of its mass flow rate (m). From this 
information enthalpies (h) may be determined, and thus the total energy flow, 
n^(hoaUet ' ^iniet)* dcHvered from the combustion Tgases, when also corfsidering 

25 incidental losses, may be determined. 
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The measurement of the energy flow of the input fuel requires knowledge 
of the heating value of the fuel and its mass flow rate. For thermal systems using 
gaseous fuels, the fuel's composition may be well characterized, thus its heating 
value may be determined based on known heats of combustion associated with 
S individual components. 

However, there are numerous situations where a fossil-fueled system's fuel 
energy flow is not well characterized. For example, ^ven a gas- fired system, 
having no on-site fuel gas analyzer, may receive fiiel from multiple sources whose 
composite heating value variation is ±5 percent or greater. The measurement of 

10 fuel flow may often present a problem of measurement accuracy, -especially at 
lower loads where flow measurement devices are not calibrated. Ineither-case the 
determination of fuel energy flow is critical for proper thermal understanding of 
any fossil-fired system, either for direct confirmation of computed results and/or 
for improving system thermal efficiency. 

IS The importance of accurately determining pollutant <:oncentiations and 

their flow rates is also critical to the practical operation of any fossil-fired system 
due to environmental cons^ints imposed through regulation, the potential of 
regulatory induced fines, and concern by the owner of the facility for 
environmental protection. 

20 Given these considerations, it is equally important to have analytical 

evidence of the errors made in the determination of fuel flow and the errors in fuel 
heating value, and thus the resultant errors in the determined thermal efficiency 
of the system. Further, any error in effluent flow, which is proportional to any 
error made in fuel flow, is significant when monitoring and reporting to regulatory 

25 agencies the effluents from any source of fossilxombustion. The method of this 
invention provides a technique for specifying and correcting pdssible^rrors biased 
on the consistency of the primary measurements of-e'ffluent C02» indicated 
Air/Fuel ratio, assumed or measured effluent HjO, possible air pre-heater leakajge 
and the t^oncentration of in combustion air. 

30 The measurement of fuel flow has traditionally been accomplished via 

measurement of its mechanical effects on a device, ^uch effects include the 
pressure drop across nozzles or orifice plates, unique fluid densities, integrated 
weighing of a fuel handling conveyor belt (conmioniy used for bulk fuels such as 
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coal), speed of sound, nuclear resonance, change in bulk storage levels, etc. 
Present industrial techniques for measuring gas or oil fuel flow result in typical 
errors of 1 percent to 10 percent relative to true values, depending on the care 
taken in designing, manufacturing, installing and calibrating the flow metering 
5 equipment and in its data reduction. Under ideal circumstances, tighter accuracies 
(i.e., smaller errors) are possible for gas and oil fuels, reaching at best 0.25 
percent, but this is considered very unusual, always requiring extraordinary 
expense. 

For bulk fuel such as coal, bio-mass, slurry fuels, wood, agricultural 

10 byproducts such as shells from nuts, trash and refuse, the typical accuracies of 
flow metering range upwardly firom 5 percent and higher. Historically, bulk flow 
measurements have such poor accuracy that they are used only as a relative 
indicator of fuel flow. For fossil-fired systems any fuel flow error greater than 
approximately I percent, and certainly greater than 2 percent, is sufficiently high 

15 to preclude trending of the monitored fuel How rate for reasons of thermal 
efficiency or for detecting d^;raded equipment. Improvement of efficiencies in 
a thermal system is classically concerned with a number of small incremental 
improvements, ^pically -each in the range of 0.2 percent to 0.6 percent. A dozen 
or more of these, taken together, may result in 3 percent to 6 percent 

20 improvement. For example, an average 4.5 percent improvement has been 
physically demonstrated at over two dozen conventional power plants, see F.D. 
Lang, "Methodology for Testing and Evaluating Power Plants Using Computer 
Simulators", 1990 Performance Software User's Group Meeting, May 1-4, 1990, 
St. Louis, sponsored by EI International, now Scientech Inc. of Idaho Falls, ID. 

25 Prior approaches which attempt to address the accurate determination of fuel flow 
are discussed below. 

Another critical consideration in determining thermal efficiency is the 
variation in the fuel's heating value due to variations in fuel chemistry. Chemical 
variations appear through the mix of fuel water, fuel mineral matter ^called fuel 

30 ash), and the relationships of the elements comprising the basic hydrocarbon 
molecular chain and any free inorganic elements: nitrogen, oxygen, carbon, 
hydrogen and sulfur, but principally carbon, hydrogen and oxygen. If an accurate 
and direct flow measurement of bulk fuels is not practical, the only alternative is 
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the determination of fuel energy flow, which is the product of flow rate and 
heating value of the fuel, based in part on the nieasurcd energy flow to the 
woricing fluid. If errors exist in the heating value, either an assumed, measured, 
or calculated value, errors will then result in the fuel flow. Prior approaches 
5 which attempt to address the determination of fuel energy flow are discussed 
below. Further, over the past 140 years of producing safe high pressure steam for 
society, mostly from coal, there has been no invention or process even -suggesting 
an ability to determine a coal's fuel ash content in r«al time based on 
thermodynamics. 

1 0 Hie approach of this invention is a much improved "Input/Loss^Method". 

Prior input/loss methods have been known to the inventor, and to T. Buna as •early 
as 1955. 

One prior approach related to the present invention was developed by T. 
Buna in 1955 for the analysis of multiple fuels fed to a power plant. His approach 

15 was to characterize a fuel's effluent COj^^given differing effluent O2 values, by 
assuming fuel chemistry of the individual fuels. With this data for multiple fiiels 
and knowledge of the Useful Energy Flow Delivered, he advocated determining 
each fuel's flow rate. He presents an "output-loss" and an "input-loss" approaches 
to determining boiler efTtciency. The present invention is related only in that a 

20 course rev»:%d from Bima's method is accomplished. This invention computes 
fuel chemistry based on effluent measurements, it assumes that all -secondary 
fuels, unlike Buna, are known having defined chemistries, heating values and As- 
Fired flows. The reference is: T. Buna, "Combustion Calculations for Multiple 
Fuels", ASME Diamond Jubilee Annual Meeting, Chicago, ILj November 13-18, 

25 1955, Paper 55-A-185. 

Another related art to the present invention was developed by the Electric 
Power Research Institute (EPRI) at the Morgantown power plant, a coal-lired 
conventional system. This technique, termed the "Output/Loss" Method, is 
described by E. Levy, N. Saninac, H.G. Grim, R. Leyse and J. Lamont, 

30 "Output/Loss: A New Method for Measuring Unit Heat Rate", Am. Society of 
Mech. Engrs., 87-JPGC-Pwr-39. The Output/Loss Method produces boiler 
efficiency ("Hijoap,) independent of fuel flow. Assuming a conventional power plant, 
by determining the energy flow to the working fluid (£mAh) and^ross-electrical 
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power production (P), system thermal efficiency may be determined, i.«., 
HboncT^/C^n^Ah). In addition, although it is not the objective of the Output/Loss 
Method, if heating value and boiler energy credits (HHV? + HBC) are known, 
then the As-Fired fuel flow (m^p) may be determined as m^p = 

5 (EmAh)/[T}boiicXHHVP + HBC)]. Useof boiler energy credits (HBC), and gross 
power (P) versus net power, are discussed in the Preferred Embodiment. The 
technique relies on measuring emission gas flow directly, and thus ilboBer 
Knowing emission gas flow allows the detennination of the majority of the 
thermal losses associated with combustion, -called ''Stack Losses". However, this 

0 approach has drawbacks when it is applied for practical applications in power 
plants for the following reasons: 1) it does not rely on measurement of flue gas 
concentrations (thus changes in the in-flows of water/steam to the combustion 
process, or changes in effluent CO2, as might reflect changes in fuel chemistry or 
air pre-heater leakages); 2) the errors in effluent -gas flow measurements in 

S irregular ducts not designed for accurate flow measurements, which is the case at 
most power plant facilities, may easily exceed ±20 percent, resulting in over ±4 
percent error in system efficiency since typically stack losses are -ZOpercent <i.e., 
±4 percent error in fuel flow); 3) the technique of direct flue ^as flow 
measurements does not consistently meet current U.S. Environmental Protection 

0 Agency's accuracy requirements of ±15 percent; 4) it is obvious Trom the above 
discussion that if the fuel's heating value (HHV) is variable, as is common with 
most coal-fired systems, and is not properly monitored in a continuous manner, 
then calculated fuel flow will also be in error, due to errors made in the assumed 
heating value; 5) direct measurement of effluent flow commonly involves 

5 ultrasonic, sonic or direct gas velocity measurements, requiring assumptions as to 
effluent compositions, i.e., fuel constituents and air pre-heater leakage, errors in 
these will force errors in the deduced efiQuent flow; and 6) if the fuel bears highly 
variable quantities of ash (a pure dilutive or concentrative influence on fuel 
heating value, but affects effluent flows through Air/Fuel relationships), then the 

0 computed fuel flow, since it is dependent on an assumed heating value with an 
assumed ash content, ^11 likely be in error 

Another related art to the present invention was presented at a technical 
conference in 1988 by S.S. Munukutla, et al. In this work the authors develop a 
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process which used effluent measurements to determine the Moisture-As-Free 
(MAF) composition of coal. Their published work teaches little relative to the art 
of monitoring thermal performance and the determination of coal chemistry of 
practical operating power plants. Munukutia, et al, do not consider air pre-heater 
S leakage, but measure combustion effluents before the air pne-heater (at the 
economizer outlet). They invoke high accuracy effluent measurements afTorded 
with gas chromatography, not common industrial instruments. This basic 
technique, commonly referred to as Thermal Analysis using in part ^as 
chromatography, is used in laboratories to detennine heating values under highly 

10 controlled conditions, refer to Chapter 9, "Measurement of Heat of Combustion" 
contained in Steam. Its Generation and Use . 40th Edition, edited by S.C. Stultz 
and J.B. Kitto, published 1992 by the Babcock & Wilcox Company, Barberton, 
OH. It was also the subject of inventor's earlier US Patent 5327356. They do not 
consider the injection of limestone as is common with fluidized bed combustors, 

1 5 representing the largest single type of steam generator sold outside the U.S . Hiey 
make no claims for the on-line determination of fitel ash, but assume it is constant 
based on laboratory analysis. Enors are made in their Equation 3, which if 
applied to their process would yield ridiculous results. This error is corrected in 
subsequent works by Munukutla as presented in 1989 and later. The ratio of 

20 atmospheric nonOj (principally N2 and Ar) to Oj is set as a constant in their work 
at 3.76 (which assumes Oj at 21.0 percent). First, this ratio, commonly found in 
fundamental thermodynamic text books, is in error; the correct ratio is 3.7737, 
herein termed ^^^^ derived from a value of 20.948 percent atmospheric Oj per 
NASA (reference U.S. Standard Atmosphere 1976, NOAA-S/T-76-1362-NAS A). 

25 Not using a correct value results in an error in determining excess air and molar 
oxygen balances. Second, in approximately 10% of the power plants tested by the 
inventor, atmospheric O2 may be degraded, resulting in a higher 4>Act value, due 
to in-leakage of flue gases into the combustion air stream and/or local atmospheric 
inversions which may result in starving the local environment for oxygen. In the 

30 present invention this ratio is expressed as a variably <|)ac» to be set by the user 
based on circumstances local to the thermal system. 

At a systems understanding level, the 1988 approach by Munukutla, et aly 
is lacking for they compute boiler efficiency only as a final result of the technique 
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using the Input*Output Method (discussed below); see their Equation 10. A 
computed boiler efificiency incorporated within the technique proper is lacking, 
thus integral consistency required for, and a feature of, the present invention is not 
assured using their methods. A most important attribute of the present invention 
5 demonstrates an integral relationship between -coal 'constituents and the 
determination of boiler efficiency, leading to system efficiency; indeed without 
coal constituents boiler efficiency simply cannot be computed by the present 
invention. Munukutla, et al., technique determines the total fiiel enei^ flow from 
system energy flow balances (their Fig. 2 and Equation 4), developing a "fuel 

10 energy" term, Am^AF " Qrad Q«eam; offering no explanation as to how the term 
"A" in this expression is related to coal heating value, as one might assume. 
However, this expression is then used to develop mass flow rates of the coal's 
constituents, termed jij (in pounds/second). These terms are then used in a 
correlation which develops total fuel energy flow (Btu/hr), thus bypassing a boiler 

IS efficiency calculation. Air pre-heater leakage is not considered in their system 
energy balances. It is noteworthy that this eneigy balance technique, leading to 
their Equation 4, is not referenced in later publications on the same subject by the 
same principal author (in 1989, 1991 and 1995), however the basic 
stoichiometrics, corrected, are referenced. The reference for this work is S.S. 

20 Munukutla, G. Tsatsaronis, Y. Shih, D.E. Anderson and S.M. Wilson, "A 
Microcomputer Software for On-Line Evaluation of Heat Rate", ASNffi Power 
Generation Conference, Philadelphia, PA, September 25-29, 1988, Paper 88- 
JPGC/PTC-1. 

In three later works, Munukutla and his colleagues attempt to refine the 
25 technique. At a 1989 conference, published in 1990, Munukutla and A. Bose 
again propose using gas chromatography to measure economizer outlet gaseous 
effluents in determining MAF heating values. The heating value correlation 
employed, using weight fractions of coal's ^constituents, was the Mott-^pooner 
correlation based on Dulong's formula. These are both well known correlations 
30 in the industry, howevor they are not based on chemical binding eneigies as this 
invention employs. However in this 1989 work, the correlation was used 
incorrectly; the correlation was intended for MAF fuel -constituents, Munukutia 
and Bose inconsistendy used the wet As-Fired. As witness to the inaccuracy of 



wo 99/49399 



PCT/US99/0426S 



-8- 

their work, their calculated MAF heating value reasonably agreed with the As- 
Fired (wet). In 1991 , Munukutla, P. Chodavarapu and D.C. O'Connor published 
essentially the same work as in 1989» again based on^as chromatography effluent 
measurements. In their 1991 work fuel water is either measured or deteraiined by 
5 difference, assuming the remaining effluents are measured on a wet base. It is 
obvious &om their work that reasonably accurate molecular weights are not used; 
given the sensitivity of any such approach, such 'errors amount to «!% enror 
(using their nummcal examples, the fuel carbon fraction is computed as 0.6616, 
versus the correct 0.6570 as determined by simply ming correct molecular 

10 weights. Further, as with the 1988 and 1989 work, no provisions are made for 
limestone injection. Again the 1991 work used the Mott-Spooner -correlation 
based on Dulong's formula to determine the heating value, HHV^^p^ but this time 
the authors used the consistent MAF fuel constituents, later correcting HHV,^p 
for fuel water and ash. The result was in error with the reported value by 387 

15 ABtu/lbm, or 3.1%. Using only their 1991 repcMed fiiel analysis, this invention's 
methods produced an error of 162 ABtu/lbm. They also reported results 6f 14 
tests preformed, resulting in a standard deviation of :<:78S ABtu/lbm, or ±5:6%. 
In 1 995, Munukuda and F. Khodabakhsh published asimilar work as in 1 989 and 
1991, but this time based effluent measurements on continuous emission 

20 monitoring system (CEMS) instrumentation. CEMS instrumentation is required 
by the US EPA on all stationary sources of pollutants. Effluent measurements 
include COj, CO, Oj, SO^ and volumetric flow rate. Under US and European 
regulations these measurements must be made at the boimdary of the system, at 
the air pre-heater's outlet, i.e., the smoke stack, not at the boiler's outlet (i;e., the 

25 air pre-heater inlet); refer to FIG. 6. Munukutla and Khodabakhsh applied the 
same stoichiometrics for this work was used since 1988, applicable for a system 
without air pre-heater leakage. Given the measurements are made with air leakage 
present, the method is fundamentally flawed; all effluents emanating firom tibe 
combustion process are of course diluted with air leakage. However, the^stem 

30 was defmed in such a manner that air leakage was factored out when considering 
a system energy balances. This was done by using the air pre-heater's hot outlet 
flow to the boiler, per their Figure 2, not the combustion air inlet to the air pre- 
heater. The consequence of this is that although fiiel flow could be computed 
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correctly, the integral determination of boiler efficiency, as with the 1988 woric, 
■ is not possible. Further, since the proposed effluent measurements in this 1995 
work are inconsistent with its basic analytics, the computed tuel constituaits are 
then seriously flawed. Limestone injection is not considered. Fuel ash is assumed 
5 constant. In all of these works, from 1988 through 1995, no mention is made of 
converting ttie resultant heating value from a constant volume base to a constant 
pressure base as required given that laboratory determined coal heating values 
derive from "bomb calorimeter" devices, a constant volume process. However, 
arefmement described in tiieir 1995 woiic suggests an iterative correction of the 
10 measured CEMS volumetric flow witii that computed through stoichiometrics. 
The present invention does not employ effluent flow measurements and thus this 
reflnement is not material. 

References for the 1989, 1991 and 1995 works by Munukutia and his 
colleagues include: S. Munukutia and A. Bose, "On-Line Elemental Analysis of 
15 Coal Using Gas Chromatography", 1989 €PRI Heat Rate Improvement 
Conference, May 3-5. KnoxviUe, TN. published 1990); S. Munukutia, P. 
Chodavarapu and D.C. O'Connor, "On-Line Coal Analysis from Measurement of 
Flue Gas Components", ASME International Power Generation Conference, "San 
Diego, CA. October 6-10. 1991, Paper 91-JPGC-Pwr-17; and S. Munukutia and 
20 F. Khodabakhsh, "Enhancement of Boiler Performance Evaluation Methods Using 
CQ^S Data", ASME International Joint Power Generation Conference, 
Minneapolis, MN, Oct. 8-12, 1995, PWR-Volume 29, 1995. 

Another approach was developed by Uie present inventor and was 
described in United States Patent 5,367,470. The inwntion of the '470 patent is 
25 noteworthy because it addresses the d^imination of boiler efficiency without 
knowledge of fuel flow and without knowledge of effluent flow, but knowing 
principally effluent CO,, effluent HjO and routine system data, and it is useful in 
many situations. However, Uie Air/Fuel ratio is not required in ttie metiiod of the 
'470 patent A key to tiie invention of the '470 patent is its requirement to 
30 repetitiously adjust, or iterate, on "an assumed water concentration in the fuel until 
consistency is obtained between the measured CO^ and Hp effluents and tiiose 
determined by stoichiometrics based on the chemical concentration of die fuel". 
Some aspects of die invention are dependent upon high accuracy effluent water 
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and carbon dioxide concentration measurements, or "predetermined accuracies" 
of these effluents. The present invention has no such limitations on accuracy. The 
difficulty in the method of the '470 patent lies with the fact that adjusting fuel 
water, which will of course alter the computed effluent water, has no prima facie 
5 effect on a dry-base effluent COj. It is true, for example, that if fuel water is 
increased, the relative fraction of the other fuel's constituents, per unity mole of 
total As-Fired fuel, will decrease assuming that the fuel's other constituents, 
nitrogen, oxygen, carbon, hydrogen, sulfur and ash, remain proportionally 
constant to each other. However, it would be unusual that any given tiiel water 

10 adjustment would produce an exactly consistent effluent COj, with the exception 
where the dry chemistiy is constant or highly predictable. Further, if the fuel has 
a variable ash content, ash having a pure dilutive or concentrative influence on 
fuel chemistry and fiiel heating value, then^uch variable effect could not possibly 
be determined by merely iterating on fuel water. A higher assumed fuel water 

IS may decrease a wet-base effluent CO2, but the actual fuel could contain much 
lower ash, thus actually increasing the amount of fuel carbon relative to the whole. 
The approach of simple water iterations of the '470 patent is useful in many 
situations, such as where the coal fuel bears little and^onstant ash, and, tiirther, 
where very high accuracy and consistent«f¥luent CO2 and HjO measurements are 

20 made, but has limitations in other applications. 

Yet another approach was developed by the present inventor and was 
described in United States Patent 5,790,420. The invention of the '420 patent is 
noteworthy as it extends the approach of the '470 patent to include -combustion 
turbine systems. The '420 patent is concerned with methods for improving 

25 thermal efficiency, determining effluent flows and determining fuel flow of fo&il- 
fired systems through an imderstanding of the total energy flow (fuel flow rate x 
heating value), the As-Fired input to the system. The '420 patent states that 
method errors will offset: "the sign of the error introduced by the heating value 
will always have an opposite change in thexalculated fuel flow". Errors may be 

30 introduced in the '420 patent by the use of its Equations 31 or 32 to<:ompute the 
dry-base heating value, dependent on knowledge of the dry molar composition of 
the fuel. These dry compositions may be determined through "use of a correlation 
relating carbon, hydrogen, oxygen and sulfur contents to a dry-base heating value 



\ 
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then correcting for water". Fuel ash is suggested by the *470 patent and the '420 
patent as being treated as a constant value. This is the case since the fuel 
constituents are solely defined per unity moles of dry-base fuel or as an As-Fired 
(for wet-base) fuel. A dry-base fuel contains ash; a wet-base fuelxontains ash and 
5 water. If the effects of variable ash were to be addressed, fuel constituents would 
by necessity initiate from a base free of both ash and water, i^e., so-called 
"Moisture-Ash-Free"; such a base is not mentioned in either the '470 patent or the 
M20 patent. The '420 patent explains that the molar quantity of fuel water "is 
iterated until convergence is achieved", resulting in an As-'Fired heating value. 

10 Again, as water is altered, the aggregate of all other fuel constituents are altered 
in opposite fashion to maintain a normalized unity moles of fuel. As with the 
approach of the *470 patent, the '420 patent requires high accuracy 
instrumentation, stating "the apparatus necessary for practicing the present inven- 
tion includes utilization of any measurement device which may determine the 

1 5 effluent concentrations of HjO and CO3 to high accuracy". The approadies of the 
'420 patent and the '470 patent, which are rudimentary Input^oss methods, are 
dependent on thermodynamically understanding a fossil-fired system without 
direct measurement of fuel or effluent flows, but there is room for improvement 
Another related art to the present invention involves using fuel water and 

20 fiiel ash mstruments intended for on-line operatioa By fuel ash is meant the fuel's 
non-combustible mineral content, before firing. These instruments employ a 
variety of techniques. Fuel water instruments include: capacitance techniques, 
microwave techniques, ultrasonics techniques and IR spectroscopy. I^uel ash 
instruments include: X-ray backscatter. X-ray fluorescence, ganmia-ray 

25 backscatter, dual energy gamma-ray transmission, gamma<-ray pair production, 
natural gamma radiation, prompt gamma neutron activation analysts, laser 
spectroscope, electron spin resonance and nuclear magnetic resonance. Knowing 
the fueFs water and ash contents can lead to adjustment of the assumed heating 
value. This provided the fuel's MAF chemistry remains constant Further, these 

30 methods do not involve techniques in which errors made in the fuel ash 
measurement are off-set by fiiel water concentrations. Further, common industrial 
accuracies are no better than ±5%, thus a ±5% error in fuel heating value. Further, 
whereas the direct determination of fuel water euid ash would aid the present 
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invention as overchecks of its calculated values, such instruments by themselves 
do not provide an integrated approach to the understanding of thermal systems. 
Clearly the computation of boiler efficiency, as integrally related to system 
parameters, is not made by simple use of these instruments. 
5 Another related art to the present invention is Thermal Analysis which 

employs laboratory techniques of differential thermal analysis (DTA) and 
differential thermal gravimetrics (DTG), also termed thermogravimetry, combined 
with gas chromatography or other high accuracy gas analyzer, to determine the 
elementary analysis and proximate analysis of coal. This technique is intended for 

10 a laboratory environment, since it analyzes only gram amounts of coal. Obviously 
using tons of coal/hour by a thermal system such as a 100 Mwe power plant, 
developing a representative gram-size sample would present an insurmountable 
problem. One of the largest companies which supplies such equipment is Mettler 
Toledo, Hightstowns, NJ which has been selling such equipment since at least 

15 1980. This basic technique is mentioned in Chapter 9, "Measui:ement of Heat of 
Combustion" contained in Steam, Its Generation and Use , sited above. 

Another ^proach was developed by the American Society of Mechanical 
Engineers (ASME) and published as its Power or Performance Test Codes (PTC). 
Several of these codes discuss two methods which are relevant: the Input-Output 

20 Method invoked in PTC 4. 1 (Steam Generators), 4.4XGas Turbine Heat Recovery 
Steam Generators), and PTC 22 (Gas Turbines); and the Heat-Loss Method 
invoked in PTC 4.1 and PTC 4.4. The Input-Output Method relies on the direct 
measurement of fuel flow. For coal-fired plants or any bulk-fuel systems, this has 
no applicability for the improvement of thermal efficiency for the above-discussed 

2S reasons of inaccuracy in fuel flow measurements. The Heat-Loss Method was 
intended to address the issue of inaccurate coal flow measurenients by 
determination of unique stack and non-stack losses, thiis: fWystem 1-0 - 
£(System Losses). The difTiculty with this method, as in the case of the 
Output/Loss method, lies in the need to make two critical determinations: 1) either 

30 measurement of gaseous efHuent flow directly or accurately knowuig fuel 
chemistry leading to an effluent-tofuel flow determination; and 2) when applied 
to coal-fired systems, measuring -effluent ash flow and its associated unbumed 
carbon in the ash. As discussed, the required accuracy is not possible on a 
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consistent basis when measuring effluent flows. Further, the traditional use of the 
Heat-Loss Method requires an "accurate, simple, and ultimate analysis of the fuel 
being fired/' The difficulty in determining both "bottom" and "fly" ash flows, 
defined in PTC 4. 1 , is evidenced by the fact the procedure is not preformed on any 
S routine basis by any known coal-fired power plant. 

In summary, the approaches of the Output/Loss Method, Munukutla and 
his colleagues, the '470 patent, the '420 patent, fuel water and fuel ash 
instrumentation, laboratory Thermal Analysis, the Input-Output Method, and the 
Heat-Loss Method all have significant limitations. The Output/Loss Method and 

1 0 both ASME Methods are flawed conceptually for at least typical large coal-tired 
systems. None of these methods are applicable for on-line monitoring. The 
methods of the *470 patent and the *420 patent do not consider: 1) the 
determination of fiiel ash (intrinsically assuming a constant or known relationship 
between fuel carbon and fuel ash); and 2) the complexities of the non-water 

IS constituents, iterating simply on fuel waler. That is, the methods of the *470 
patent and the ^420 patent simply alter fuel water under the assumption that the 
relationships between the fuel's non-waterxonstituents remain as constants or are 
predictable through correlations. If the determination of the fuel's constituents is 
flawed, then the determination of the fuel's heating value is flawed, and thus the 

20 fuel flow will be in error. Although as the ^470 patent and the '420 patent^tate, 
such errors in heating value and fiiel flow tend to off-set one another, this is 
typically true only if the ash is both a relatively small fraction of the total fuel and 
of essentially constant concentration. As will be seen in relation to the present 
invention, tolerating such errors voids useful information associated with absolute 

25 knowledge of heating value, voids the accurate determination of effluent flow, and 
voids any computational ovoxheck of the accuracies of the effluent 
measurements. None of these methods, except the '470 Patent and the '420 
patent, consider multiple fuels which are commonly used in commercial power 
plants. None of these methods considers the thermodynamic determination of 

30 fiiel ash. None of these methodsconsiders the use of limestone injected into the 
combustion process. None of these methods, excqst the '470 Patent and the '420 
patent, consider air pre-heater leakage such that gas concentrations on either side 
of the air pre-heater are useable in stoichiometric relationships. None of these 
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methods, except the '470 Patent and the '420 patent, consider variable O2 in the 
combustion air local to the system. 

Another approach was developed by the United States Environmental 
Protection Agency (EPA) as related to the determinations of effluent flow and 
5 individual emission rates (lbp^,|^„/million-BtUfi^,). EPA's sqsproach is described 
in its regulations. Chapter 40 of the Code of Federal Regulations (40 CFR). 
Speciflcs are described in 40 CFR Part 60, Appendix A, Niethods 1, 2 and 17 
defining various techniques for measuring effluent flows, and in 40 CFR Part^O, 
Appendix A, Method 19 defming "F Factors" used to determine^emtssion rates. 

1 0 The EPA requires the direct measurement of effluent flow from stationary sources 
of fossil combustion. The EPA also requires the reporting of emission rates for 
the major pollutants, determined using the F Factor technique. The EPA's 
approach has the same shortcomings as discussed for the Output/Loss Method. 
The EPA approach does not require any inter-relationship between a x:omputed 

15 fuel flow, which through stoichiometrics must be consistent with effluent flow. 
In a work by Lang, et aly reporting actual test results on a large power plant, 
EPA's Methods 1, 2 and 17 resulted in -effluent flows, on average, 12 percent 
higher than those consistent with system efficiency; this implies a 12% higher fuel 
flow, a large error in system understanding. Further, all EPA methods produced 

20 higher flows relative to those consistent with system efficiency, and none were 
consistent in themselves. These results are typical of the 12 to 13% bias reported 
throughout the industry. See F.D. Lang,^/ al, "Confirmatory Testing of the 
Emissions Spectral Radiometer/Fuel Flow (ESR/FF) Instrument", Electric Power 
Research Institute (EPRJ) 1994 Heat Rate Improvement Conference, May 3-S, 

2S Baltimore, MD. In another study at two large power plants, use of EPA's Method 
2 produced 9.8 percent and 18.6 percent higher system heat n^s (i.e., system 
efficiencies) based on measured effluent flows, see R.D. McRanie, et cdy "The 
Electric Power Research Institute Continuous Emissions Monitoring Heat R^^e 
Discrepancy Project, An Update Report - December 1996", available from EPRI, 

30 Palo Alto, CA. Studies by Lang and M.A. Bushey, analyzing 14 power plant test 
results, indicated that errors in emission rates using the EPA F Factor method 
could range from -8 percent to +4 percent; and, when studying Tive dozen coal 
samples, they found that for 18 percent of the samples the error exceeded 5 
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percent. See F.D. Lang and M.A. Bushey, "The Role of Valid Emission Rate 
Methods in Enforcement of the Clean Air Act", EPRI 1994 Heat Rate 
Improvement Conference, May 3-5, Baltimore, MD. 

Another related art to the present invention was developed by Roughton 
5 in 1980; see JJE. Roughton, "A Proposed On-Line Efficiency Method for 
Pulverized-Coal-Fired Boilers", Journal of the Institute of Eneigy, VoK20, March 
1980. pages 20-24. His work served in part as the bases for the abovesited EPA 
methods, and is related to the Output/Loss Method Roughton's method prodiiccs 
boiler efficiency {r\^^ independent of fuel or effluent flows. His work computes 
10 boiler efficiency from the process* stack and non-stack thermal losses, ^tvaluated 
per imity of As-Fired fuel flow. Of these losses, the major loss is the stack loss. 
Whereas this is directly measured for the Output/Loss Method, Roughton relies 
on the statistical relationship between dry effluent flow and total As-Fired fuel 
energy flow assuming a water-free (dried) fuel. He relies on an observed ratio of 
1 5 these two quantities being essentially constant at 0.0008257 Ibm^yen/Btu^^j^ Fod 
(referred to below by the term Lf^, and called the "fuel factor^O- The method's 
accuracy is based solely on this value remaining constant The EPA's F-Factor 
approach fundamentally relies on this same ratio remaining constant, see the work 
by Lang and Bushey. It has been found that for a specific fuel, having a^ertain 
20 Rank from a certain geographical region, this ratio is indeed constant; but found 
not the same for different fuels. Further, it makes no claim as to predicting 
heating value, indeed Roughton states: "Using this ratio it is possible to obtain the 
dry stack loss and moisUire loss without the need for an ultimate analysis of the 
coal or for an accurate determination of calorific value" of the As-Fired fuel. 
25 Complete thermodynamic understanding of fossil-fired systems, for the 

purposes of improving system efficiency and regulatory reporting, requires the 
determination of fuel mass flow rate, fuel chemistry, fuel heating value, total 
effluent flow rate, emission rates of the common pollutants, and thermal efficiency 
of the overall combustion process. All such quantities must be <ietOTnined with 
30 thermodynamic consistency. There is a need for an improved approach to these 
detenninations. 
/// 
/// 



W0 99/493W 



PCT/liS99/04265 



-16- 

SUMMARY OF THE INVENTIOlvJ 

The approach of the present invention, termed the Input/Loss Method, 
consistently detennines fuel flows, effluent flows, emission rates, fuel chemistry, 
fuel heating value and thermal efficiency, resulting in improved determinations 
5 of the thermal efficiency of any fossil fueled system. The Input/Loss Method has 
been applied through computer software, installable on a personal computer, and 
demonstrated being fully operational. This computer and software is termed a 
Calculational Engine, receiving data fipom a fossil fueled system's data acquisition 
devices. The Calculational Engine operates continuously, i.e., in "real time" or 

10 "on-line", as long as the fossil fueled system is receiving and burning fuel. 

Prior to on-line operation, the Input/Loss Method requires certain 
initializing data involving reference fuel chemistry and heating value, and 
reference fuel stoichiometric data associated with the reference fuel chemistry. 
In addition, those computer programs which will describe the steam generator's 

1 S air handling equipmrat and its heat exchangers, or the ^ turbine and its heat 
exchangers and/or the steam turbine cycle, require routine initiating input. If 
operating for the first time, initial estimates of fiiel chemistiy and heating value 
are required, typically the reference values are used. Effluent measurements are 
required If off-line, the assumed effluent measurements may be either consistent 

20 with the assumed fuel, or take a bias for the study of instrumentation^rror. Using 
these data, error analyses are preformed indicating which assumptions will yield 
minimum errors by exercising the Input/Loss Method as an analytical tool. Two 
such assimiptions are required. The first of these is how the molar fraction of fuel 
ash, using a Moisture-Ash-Free (MAF) base, aMAF-Asip should be treated: using a 

25 constant value of Omaf-aji^ ot correlating awAF-Ash as a function of MAF heating 
value, HHVmaf; measuring it directly using a fuel ash meter, or, preferably, 
determming a^AF-Ash by explicit solution requiring the measurement of the 
system's wet combustion Air/Fuel mass ratio (a relative measurement routinely 
made in all fossil fueled systems). In general, this Air/Fuel ratio requires 

30 normalization given measurement bias such that stoichiometric consistency is 
achieved. An alternative to the Air/Fuel ratio is the indicated fuel mass flow rate, 
also routinely measured in all fossil fueled systems, but requiring normalization 
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(using error analysis procedures of this invenUon). The second assumption is how 
the molar fraction of fuel water, using a MAF base, o^.^ should be inOuenced 
by effluent water: using a constant value of the effluent H|0. is preferred; or the 
effluent stream may be instrumented for the direct measurement of HjO. 
5 When operating in real time and using the initialized data, the Input/Loss 

Method performs the following sequential steps: 1 ) obtain measuienients of the 
concentrations of the common pollutants to accuracies^nimon to the electric 
power industry; 2) obtain measurements of the gross shaft -electrical power, 
mechanical power, and/or the useful energy flow developed firom the system; 3) 
10 if multiple fossil fuels are used, their properties are combined {e.g., using the 
FUEL.EXE program) to form a composite fuel, composite higher heating value, 
and, even if a single fuel, to prepare input for step 4; 4) fuel concentrations and 
heating value are input to a steam generator or gas turbine computer simulator 
(e.g., the EX-FOSS£XE program); 5) obtain measurements of the effluents O,, 
15 coj, SOi. and H^O if appropriate, and the indicated Air/Fuel ratio (or indicated 
fuel mass flow) if appropriate, to accuracies common to the electric power 
industry: 6) the steam generator, or^as turbine, computer simulator is executed 
producing consistent stoichiometrics given the supplied fuel and input of tiie 
measured effluent O, and common pollutants, including the computed effluents 
20 COj, SOj and HjO values, the Air/Fuel ratio, the moles of fuel per 100 moles of 
dry gaseous effluent (termed x). and at least the following thermal performance 
parameters: fuel and effluent floM«. system thermal efficiencies and emission rales 
- all consistent with the input fiiel's chemistry and heating value; 7) by solution 
the molar MAF fractions of fuelcarbon, water and sulfur arecomputed as explicit 
25 stoichiometric solutions, both dependent principally on the effluents O^, C(^,^02 
and HjO (which was not done in the '470 patent or the ^420 patent); 8) through 
dependency on the molar MAF fraction of fuel carbon, the molari^AF fractions 
of fuel nitrogen, oxygen and hydrogen are tlctermined; 9) as optioned in the 
hiitialization, the molar MAF fraction of fuel ash is determined, the preference is 
30 by explicit solution; 10) the molar MAF fuel species are converted to a molar dry 
base, titen converted to a molar wet base (As-Fired), and finaUy to wet weight 
fractions (As-Fired); 1 1) the higher heating value is -computed based on changes 
to the fuel's MAF constituents, tiien converted to a dry base, and then to an As- 
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Fired base; 12) the resiilts of the last two steps, fuel chemistry and heating value, 
are then input to the FUEL.EXE program (or a similar program) of step 3 and the 
processes repeated until convergences on the fuel moles (x), As-Fired heating 
value (HHV^f)» f^^l flow and minor stoichiometric terms are achieved; 1 3) after 
S convergences are achieved, and if errors are within criteria, the procediure is 
deemed successful, if not instrumentation is identified by the process allowing for 
both correction and the minimization of errors through application of 
multidimensional optimization techniques; and fuially, 1 4) after convergences are 
achieved and error analysis completed, the operation of the system is adjusted to 

10 improve its thermal efficiency and/or to minimize die polluting emissions and/or 
to report effluent flow and emission rates to regulatory authorities. When this 
process is completed all objectives of this invention will have been met. 

The fossil fueled system operator has assurance of complete 
themiodynamic understanding of the system because of: 1 ) explicit relationships 

IS between measured effluents and the key fuel constituents of carbon, water and 
ash; 2) an explicit relationship between these and the computed heating value; 
3) an explicit relationship between the fuel eneigy flow (heating value and flow) 
and the useful energy flow developed from the combustion ^gases; and 4) an 
explicit relationship between fuel flow and effluent flow. 

20 The apparatus necessary for practicing the present invention includes any 

measurement device (or combination of devices) which determines the effluent 
concentrations of Oj and COj,, and, if appropriate, effluent H2O, CO & SOj to 
current standards found in the electric power industry. Further, the system's 
routinely measured Air/Fuel ratio (or fuel mass flow) is required provided it is 

25 consistent, not necessarily accurate (normalization of the signal is provided) if fuel 
ash is to be determined in real time. Further, routine and common thermal system 
data, all of which are typically known to thermal system operators, is required 
such as: effluent gaseous temperature; combustion air psychrometrics; working 
fluid flows, pressures and temperatures at key heat exchangers, and the like. 

30 It is therefore an important object of the present invention to provide a 

procedure for determining the energy flow of the input fiiel to a fossil fueled 
system without direct measurement of the fuel flow rate or heating value or fuel 
chemistry, and in accomplishing this, to assure that system efficiency, and system 
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mass and energy in- flows and out-flows are consistent. 

It is a further object of the present invention to provide a quantitative 
procedure of demonstrating the consistency of the Method's results. 

It is a further object of the present invention to provide a means Tor 

5 determining the energy flow of the input fuel of a Tossii fueled system by 
predicting the composition of the input fuel, including its ash content, and with 
this information predict heating value, and then back-calculate the inputluel flow 
rate from a classical use of system efficiency. 

It is a further object of the present invention to provide a means for 

0 determining both the total effluent flow rate (cubic feet/hour or pounds/hour), the 
emission rates (pounds/million-BtUfuei) rates (pounds/hour) of all 

effluents including the common pollutants produced from a fossil Tiieied system 
by determining the fiiel flow rate indirectly and having knowledge of the fiiel's 
chemistry and effluent concentrations. 

S It is a further object of the present invention to provide a procedure for 

determining the thermal efficiency of a fossil fueled system without directly 
measuring the input fuel flow rate. 

It is a further object of the present invention to provide an intrinsic self- 
checking procedure of the Method of this invention, computed in real time, in- 

0 process, in which the computed and measured effluent CO2, -effluent HjO, 
Air/Fuel ratio and computed fuel factor are compared for-consistency, and if errors 
are within criteria the procedure was deemed successful, or not. 

It is a further object of the present invention to provide a procedure to 
identify which of the effluent CO2, effluent HjO and Air^uel ratio measurements 

S is producing erroneous measurements, such that corrective actions may be taken 
if appropriate. 

It is a further object of the present invention to provide a procedure in 
which the errors made in predicting the concentration of fuel water will be off-set 
by fuel ash, though an explicit computation of fuel ash. Thus any enror in fuel 
0 water, as off-set by fuel ash, will have negligible effect on the As-Fired heating 
value. 

It is a further object of the present invention to provide a means for 
determining the flow of a fuel's ^olid nonnrombustible mineral material 



wo 99/49399 



PCT/US99/042$5 



-20- 

(commonly referred to as fuel ash) associated with coal fuel, and thus adjust and 
improve the operations of the ash removal equipment in a combustion system's 
effluent stream. 

It is a further object of the present invention to demonstrate that all 
S aforementioned objectives associated with a fossil fueled system are also, and 
herein declared, objectives associated with any fossil fueled system producing 
electrical power, mechanical power and/or useful energy flow from the system. 

It is a further object to provide an approach which yields improved results 
as compared with EPA Methods 1, 2, 17, and/or 19. 
10 It is a further object to provide an approadi which yields improved results 

as compared with the methods of the ASME Power Test Codes 4. 1 , 4.4, and/or 22. 

Other objects and advantages of the present invention will become 
apparent when the Method and apparatus of the present invention are considered 
in conjunction with the accompanying drawings. 

15 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram illustrating the initialization procedures 
including the study of instrumentation errors; 

FIG. 2 is a block diagram of the present approach, showing the 
determination of fuel flow, system efficiency, fuel -chemistry, fuel heating value, 
20 effluent flow and emission rates; 

FIG. 3 is a plot of Moisture-Ash-Free molar fuel hydrogen versus 
Moisture- Ash-Free molar fuel carbon using typical coal data; 

FIG. 4 is a table of suggested Calculational options associated with this 
invention allowing for the analysis of any hydrocarbon fuel and^als of various 
25 Ranks; and 

FIG. S is an error plot showing laboratory-reported wet heating values 
versus those computed by the Method of this invention for axoal-flred power 
plant having wildly varying fuel chemistries and heating values. 

FIG. 6 is a schematic representation of a conventional or fluidized bed 
30 power plant i I lustrating use of stoichiometric relationships. 
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PESCRtPTION OF THE PREFERRKD FMBODIMENT 

For clarity, note that FIG. 6 depicts fossil-fired power plant thermal system 
20. In this power plant thermal system 20, a fuel feed 22 and combustion air 24 
are provided from an upstream side 26 of the heat exchangers/combustion region 
S 28. Note diat this region 28 does not include the air pre-heater 36. In addition, 
in some types of pow^ plant thermal systems 20, other materials may be injected 
into the heat exchangers/combustion region 28, such as a flow of limestone 3 1 to 
minimize effluent SO2 by chemically binding sulfur. The fuel feed 22 contains, 
in general, combustible material, water and mineral matter (called fuel ash). The 

10 fuel ash is an unbumable component that passes through the system with little 
physical change, but which is heated and cooled. 

In the heat exchangers/combustion r^ion 28, the fossil fuel 22 is burned 
with the combustion air 24 to form hot combustion products. Heat from the 
combustion products is transferred to a working fluid 30 that flows through heat 

IS exchangers 32 that are depicted as integral with the beat exchanger/combustion 
region 28. The heated working fluid 30a is used in a manner appropriate to a 
working fluid to generate a useful ou4)Ut 33 (for aconventional power plant such 
useful output may be supplied to a turbine cycle which could produce^lectrical). 
There may be a water leakage 29 into the hot combustion products, not associated 

20 with the fuel feed 22 as a result of, for example, soot blowing or tube leaks from 
the heat exchangers. 

After leaving the heat exchangers/combustion region 28 on its downstream 
side 34, the cooler combustion products commonly flow through ducts 35 which 
may contain effluent ash removal equipment, passing then to an air pre-heater 36, 

25 where a further portion of their heat eneigy is transferred to an incoming air 
stream 38, which then becomes the combustion air 24. The total air delivered to 
20 is 25. In many cases, a leakage flow 40 of air enters the flow of combustion 
products from the incoming air flow 25 as it passes through the air pre-heater 36. 
The further-cooled combustion products leave the air pre-heater 36 and p^s to a 

30 stack 42 and are then exhausted to the atmosphere. 

This general system description is applicable to a wide variety of fossil- 
flred power plants, such coal-burning power plants, oil-buming power plants^^gas- 



wo 99/49399 



PCTAJS99/04265 



-22- 

fired power plants, biomass combustors, fluidized bed combustors, a conventional 
electric power plant, a steam generator, a package boiler» a combustion turbine, 
and a combustion turbine with a heat recovery boiler combustion turbines. This 
list is not meant to be exhaustive, however, and is presented to illustrate some of 
S the areas of applicability of the present invention. 

It is desirable to quantify the operation of this power plant thermal system 
20. If it is characterized quantitatively, then the relationship of variations in the 
operating parameters to the inputs and ouQ>uts of the power plant thermal system 
20 may be understood. This understandii^ in turn, permits the operation of the 

10 power plant thermal system 20 to be optimized for thermal efficiency and 
pollution minimization. 

In this system, some quantities are readily measured, and others cannot be 
measured on-line, in real time with accuracy sufficient to quantify the operation 
of the power plant thermal system -20 to the required accuracy. For example, ^as 

15 concenurations, working fluid flow rates, and temperatures and pressures may be 
readily measured with good accuracy by conventional sensors located at a defined 
system boundary 44 and elsewhere such as in the regions 3S and 42. On die other 
hand, the As-Fired fuel chemistry of the fuel feed 22 cannot be measured 
accurately in real time if the fuel is coal or other bulk fuel. The present invention 

20 uses a modeling analysis to quantify the As-Fired fuel chemistry of the fuel feed 
22 and other parameters in a self-consistent manner, which allows the-entire 
power plant thermal system 20 to be quantitatively understood. Once the power 
plant system is quantitatively understood, its operating parameters may be varied 
to optimize its thermal performance. 

25 The present invention provides an Input/Loss Method which allows for a 

complete understanding of fossil fueled combustion systems such as power plants, 
through application of non-direct but explicit detennination of fuel and^ffluent 
flows, fuel chemistry, fiiel heating value and thermal <et!iciency, resulting in 
improved thermal efficiency. In addition the Input/Loss Method of the invention 

30 determines the effects of inaccuracies in the Method's primary measurements. 
This invention was developed primarily for<:oal-fired systems. However, this 
invention may also be applied to any fossil fueled combustion process including, 
but not limited to, gas, oil, bio-mass, slurry fuels, wood, agricultural byproducts 
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such as shells from nuts, trash and refuse. 

More particularly, this invention relates to the Input/Loss Method for 
monitoring the operation of a fossil fueled process by analyzing the composition 
of combustion effluents COj, O2, SOj and, depending on the characteristics of the 
5 fuel and combustion process, effluent HjO and/or the systems' indicated 
combustion Air to Fuel ratio. Other routine system data are also required such as 
gaseous effluent temperature. Having computed the thermal efficiency and 
consistent fuel flow rate* and having computed the fuel's chemical composition, 
and having computed the fuel's heating value, the plant's effluent flow rate and 
10 emission rates may then be determined. Further, this invention allows for the 
determination of the effects of inconsistencies in the principal measurements of 
effluent COj, effluent l^O and the systems' Air/F uel ratio (or indicated fiiel flow), 
as such inconsistencies impact the thermal efficiency of the system and its effluent 
flow rate. 

IS The subject of this invention includes any fossil-fueled system producing 

a useful and measurable output Such fo^il-fossil systems include, for example, 
conventional electric power plants; fluidized bed electric power plants; any boiler 
whose combustion gases cause the increase in energy of a measurable fluid; a 
combustion turbine producing electrical or mechanical power; and a combustion 

20 turbine producing both electrical and/or mechanical power, and whose combustion 
gases cause the increase in energy of a measurable fluid. Although the traditional 
working fluid from fossil-fired systems is steam, this invention places no 
limitations on the nature or quantity of useful output associated with fossil 
combustion. However, if fuel and effluent flows are to be determined, the 

25 invention does require that the working fluid's energy received firom the 
combustion process be measurable. This invention applies to a wide range of 
systems, from small boilers as might be used for space-heating hotels, to large 
conunercial electric power plants. Working fluid flows, such as steam flow, could 
typically range from 10,000 pounds/hour (Ibm/hr) to 10,000,000 Ibm/hr. 

30 Electrical or mechanical power could typically range from S Mega- Watts (MW) 
to 1,500 MW. 
/// 
/// 
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System Cakul^tiQHS 

The present invention allows for a complete understandii^ of fossil-fueled 
combustion systems through explicit detemiinations of fuel and effluent flows, 
fiiel chemistry, fiiel heating value and thermal efficiency through thermodynamics, 
5 not through direct measurement of fuel flow, fuel chemistry nor through ^effluent 
flow. The approach principally relies on measurements of the effluents from 
fossil-fueled systems, and other routine system data. Given the nature and choice 
of such data, it has the potential of high resolution on a <;ontinuous basis, but the 
process only anticipates common industrial accuracy standards of required 

10 instrumentation. Hie data may be input to a computer program for resolution of 
mass and energy balances associated with the system. Measurements include the 
concentration of combustion gases exiting the system ti e., smoke stack), the 
direct production of power, and the total energy deposition to the working fluid 
from the hot combustion gases. 

15 The thermal efficiency of a fossil-fueled system, termed system thermal 

efficiency is defined in the usual manner (i.e., total useful output divided by total 
inputs): 

Tlsysicm = Otoss Shftf^ Powef (CT) + Usrful En^gy F)qw Pelj v^gd 

Fuel Energy Flow + System Energy Flow Credit ( 1 ) 

20 In Equation ( 1 ), "Gross Shaft Power (CT)" refers to direct shaft power developed 
from the combustion gases as normally associated with a combustion turbine. 
Such a system may or may not employ a heat recovery steam generator (if it does 
then the "Useful Energy Flow DeliveFed** term will be non-zo^o). For a 
conventional fossil- fueled power plant the Gross Shaft Power is zero, and the 

25 Useful Energy Flow Delivered refers to the energy flow to the working fluid of 
a turbine cycle (TC). In general the terms Gross Shaft Power (CT) plus Useful 
Energy Flow Delivered define the thermal system's "useful outpuf. For a 
conventional power plant or a conventional boiler, it is necessary to separate the 
power production process, or turbine cycle, firom the combustion process. Note 

30 that the efficiency of the power production process could be thought of as the 
efficiency associated with the working fluid's process. 
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^lurbiiMKycie Usef ul Powcr Produced (TC) 

Useful Energy Flow Delivered (2) 

Therefore: 

'Hsystein "~ ^boiler ^ttiibine<ycle (3) 

S where for a conventional boiler: 

TlboHer = UscfuI Energy Flow Delivered 

Fuel Energy Flow + System Energy Flow Credits (4) 

For a combustion turbine system Tiju,btn«yc,c is unity, and the r]^^„ term is 
identically tj^j^ as defined by Eq. (1). In all further discussions the tern 

10 is used in place of r\^^^^ when referring to a combustion turbine system 
(understanding that the value of Titu,binecyde = 1 00), and the "Useful Energy Flow 
Delivered" of Eq.(4) is replaced with "Gross Shaft Power (CT) + Useful Energy 
Flow Delivered" of Eq. (1). Further, for convenience of discussion, combustion 
turbine efficiencies employ a higher heating value; conversion to a lower heating 

1 S value base, as commonly used for-combustion turbine ^systems in North America, 
is a simple conversion that is well known in the art and described in ASME PTC 
22. Outside of North America, lower heating values are used for most analyses 
of fossil fueled systems. Several authors in this field disagree with the use of the 
term "System Energy Flo\y Credits"; see the denominator of the efficiency 

20 Equation (4). Its use, or not, is not important to this invention; however, its 
consistent treatment relative to r\^^i„ is most important If energy credits are to be 
used, >^ch is the Preferred Embodiment, then the term must appear in the boiler 
efficiency calculation, and also when fuel flow is computed, i.e., Eqs.(19B) and 
(21 ). If the term is not used, then only fuel flow and heating value ^1 define fuel 

25 en^gy flow and thus iiboiier - 

The evaluation of electrical or mechanical power terms is well known. The 
Preferred Embodiment uses gross electrical or mechanical power, that is total 
power derived directly from the combustion^ases. However, net power may be 
used, that is gross power less that power consumed within the system, called 

30 "house load". Use of either of these terms is not important to this invention, only 
that gross or net terms be used consistently. The evaluation of the term Useful 
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Energy Flow Delivered is well known, and involves simply measuring of mass 
flow rate, and pressures and temperatures of the working fluid. Therefore the 
evaluation of 'ntj,rt,i,«Hg«ie is known to those skilled in the art, and may be 
accomplished with high accuracy and reasonable expense. 
5 It is evident from the above discussion that all principal terms are 

descriptive of any fossil-fiieled system producing useful eneigy flow and/or power 
(in any form). The Method of this invention is equally applicable to a combustion 
turbine as to a conventional boiler producing useful energy flow by heating a 
working fluid such as steam or liquid water, and to a<:onventional fossil-fueled 

10 power plant consisting of a conventional boiler heating liquid water to create 
steam for direct usefulness and/or for a steam cycle producing any combination 
of electrical or mechanical power, defined as useful power from a turbine cycle. 
The Method of this invention is equally applicable when using higher or lower 
heating value as a bases for the efficiency definition.. 

15 The boiler efficiency of Eq. (4) is helpful for continuous monitoring of 

thermal performance if the fiiel flow and other parameters comprising r\^i„ may 
be measured accurately. Even if so monitoring, then the directly measured fuel 
flow, fuel water, fiiel ash and/or the resultant boiler efficiency may be»compared 
to the calculated values via procedures herein disclosed. The Method of this 

20 invention computes iiboaen although Eq.(4) (or Eq.(l) for a t:ombustion turbine) is 
reformulated such that fuel flow is calcuiationally excluded. After excluding the 
fuel flow term, three major deficiencies in the knowledge of a fossil-fueled 
system's thermodynamic process remain, and especially when considering a coal- 
fired system: 1) the complexities of the combustion process (and if axjombustion 

25 turbine, the complexities of the compression and expansion of -gases); 2) the 
specification of thermal losses not directly related to the combustion process <but 
which could affect measured fuel flow); 3) the variability or uncertainty of the 
heating value; and 4) the variability or uncertainty of the fuel chemistry as it 
impacts r\^i„ and the heating value. 

30 A computer program EX-FOSS.EXE has been developed to address the 

first two of these difficulties. It is a commercially available program which has 
been in use in the power generation industry since 1985 and is available (together 
with the FUELJEXE program) from Exergetic Systems, Inc. of San Rafael, 
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California. The methodology of the EX-FOSS.EXE program is printed in detail 
in the *420 patent, see its Figures lA and IB, and their associated discussion, 
which figures and discussion are incorporated by reference. The FUEL.EXE 
program, using both known fuels (e.g., stabilizing gas, oil and/or coals) and the 
S fuel being resolved, performs simple weighting <:aIculations resulting in a 
comi>osite fuel. FUEL.EXE may accept any number of fuels, each of which is 
described by its elementary molar or weight fractions. Gaseous -components are 
input by molecular molar or weight fractions. Each descriptive fuePs Dractions 
must sum to imity. These are combined to produce a -composite fuel, at the user*s 
1 0 option, using a molar or weight base. The individual fuers heating values are also 
weighted. This composite is then supplied to EX-FOSS.E3ffi as required input. 

The EX-FOSS methodology separates the definition of boiler efficiency 
into components which, taken separately, calculationally exclude the first of the 
aforementioned problem areas. When separated, terms called combustion 
IS efficiency, t^^, and boiler absorption efficiency, t)^, are developed. Use of r\c 
allows consideration of only input and ou^ut terms which may be measured with 
high accuracy, thus eliminating resolution of combustim's thermodynamic 
complexities. T]e addresses the major loss in a fossil fired system, i^«., stack losses. 
The use of ti^ allows for the consideration of only "non-stack" losses, losses which 
20 generally have minor impact on system efficiency, thus whose understanding does 
not require high resolution. The following leans are defined, and typical units aa?e 
provided for clarity: . 

m^F s As-Fired Fuel Mass Flow Rate; Ibm^i/hr. 
HHVP = Higher Heating Value at Constant Pressure; Btu/lbmAp. 
25 HBC = Specific System Energy Credits; Btu/lbmAp. 

EF = Fuel Energy Flow (fuel flow x higher heating value 
at constant pressure); Btu/hr. 
sm^pHHVP 

CF 5 iSystem Energy Flow Credits (fuel flow x specific energy 
30 credits relative to the system); Btu/hr. 

= m^^pHBC 

HPR s Enthalpy of the Combustion Products (includes the heat of 

formation plus jCpdT at the effluent boundary); Btu/lbm^p. 
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HRX s Enthaply of the Reactants (based on the heating value, 

sensible heating and energy credits); Btu/lbniAP 
HSL s Specific Stack Losses (includes losses directly 

affecting the energy released during combustion, defined by 
5 PTC 4.1: Lc» L,„f» ^h* Liiia» Lx> L^, Leo, & ^uhc^ 

divided by m^p); Btu/IbAp. 
HNSL s Specific Non-Stack Losses (boiler or combustion turbine losses 
whose mechanisms originate from the combustion process 
or from the hot gases, and interface directly with the 
10 environment thus a direct -effect on the BBTC term; defined 

by PTC 4.1 as: Le, Lp, L^, L, & L^jci defined by PTC 4,4 as: 
Lg & Lvv; in addition, if appropriate, to the gas turbine's 
shaft bearing losses, applicable gear losses and other similar 
fKctional losses associated with delivery of useful shaft 
15 power directly by t:ombustion gases (all divided by m^p); 

Btu/lbm^p. 

ERC = Energy Released during Combustion; Btu/hr. 
= EF + CF - £(Stack Losses) 
^EF + CF-mApHSL 
20 = mAF(HPR - HRX) 

BBTC = Gross Shaft Power (CT, if a combustion turbine) 
plus Usefiil Energy Flow Delivered; Btu/hr 
= ERC - E(Non-Stack Losses) 
= ERC-m,,FHNSL 

25 With these definitions, equivalent ways to express fossil boiler efficiency include 
the following. 

Tlboiiero BBTC 

EF + CF (5) 

r\^^^ ERC -EfNon-Stack Losses) 
30 EF + CF (6) 
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r\toi}^ . fEF + CF^ - EfStack Ln^<;^<:^ - ETNon-Stack Losses^ 

EF+CF (7) 

r\^ii„ „ ERC 11 .0 . gflN^Qin^St^ck los^e^)\ 

EF + CFl ERC J <8) 

5 Eq.(8) suggests that boiler efficiency may be divided into two 'separate 
efficiencies: one descriptive of the combustion process per se (called the 
combustion efficiency), and the other descriptive of -certain non-stack losses 
(called the boiler absorption efficiency). 

The combustion efficiency definition is suggested from the Input-Output 
1 0 Method defined in PTC 4. 1 or in PTC 4.4: that is, net aiergy flow produced at the 
thermodynamic boundary divided by the total eneigy flow input -(fuel energy flow 
and system energy flow credits). The following develops r\Q on a unity tuel flow 



Tic = 1 >0 - ^^^'^t^ Losses) 
15 EF + CF (9) 

= fPF + CF> ■ Ef^tack Loss^) ■ 

EF + CF (10) 

Tic- 

EF + CF (11) 
20 Tic = [ niAF (HPR - HRX)) / (EF + CF) (12) 

s HPR - HRX 

HHVP + HBC (13) 

The boiler absorption efficiency is derived from the Heat-Loss Methods 
found in the PTC 4.1 and 4.4 although limited to non-stack-energy terms. It must 
25 be referenced to the Energy Released during Combustion term (ERC) if all losses 
(on a systems basis) are to be additive when-calculating the boiler efficiency, see 
Eq.(7): 

Tl^ = 1 0 ■ S^op-Stac R Losses), 

Tic(EF + CF) (14) 
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The quantity nc(EF + CF) defines the ERC term, see definitions and Eq.(l 1), 
thus: 

tlA= 1.0 - m^pHNSL/ERC (ISA) 
T,^=1.0, HNSL 

5 HPR-HRX (I5B) 

The Stack Losses includes the following PTC 4.1 terms relating Stack 
Losses to total As-Fired fiiel flow rate: 

mApHSL = Lq + L„F + L„ + L^ + Lx + L^ + Leo Luh + Luhc (16) 

The quantity HNSL includes die foHowir^ PTC 4. 1 terms, PTC 4.4 terms, 
10 and terms relating to turbine losses (energy flow losses in deliveringgross shaft 
power): 

mA|HNSL = La4j + Lp + Lu + L^ + Luc +Ihm.4 "^W + (turbine shaft losses) 

(17) 

The combination of the combustion efficiency and boiler absorption 
15 efficiency is equivalent to the classical definition of boiler efficiency as defined 
by ASME Power Test Codes (e.g., PTC 4.1), British Standards (e.g., BS2885), 
and others. The following, using direct energy flow tains, ^ain demonstrates the 
derivation of boiler efficiency, see Eq.(13) and Eq.(lS): 

nboUcr^^^C^A (18) 

20 nboaef= HPR-HRX . [HPR-HRX-HNSL]/(HPR-HRX) (19A) 

HHVP + HBC 

HHVP + HBC (t9B) 

Tlboi!cr = 



25 m;,p(HHVP + HBC) (20) 

Eq. (20) may be rearranged to solve for fuel flow rate, provided heating value and 
BBTC are accurately known. The resolution of fuel flow and heating values are 
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objectives of this invention: 

BBTC 

t1j„„«(HHVP + HBC) (21) 

By separating boiler efficiency into combustion and boiler absoiption 
5 components, the analyst has knowledge as to where degradations 2ire occurring. 
If combustion efficiency decreases (i.e., stack losses increase), the plant engineer 
would consider: fuel-air mbcing equipment, degradation of hardware directly 
interfaced with delivery of the combustion air, low beat content in the fiieU 
improper operation of the fuel's burner mechanism, improper operation of the 

10 combustion proper such as the location of the fire-ball, etc. - all sources directly 
affecting the combustion process (i.e., stack losses). The terms comprising 
combustion efficiency may be easily reduced to a unit basis of total as-'ficed fuel, 
refer to Eq.( 1 3); as such these tenns have the potential to be determined with great 
accuracy. HHVP is the higher heating value; corrected for an assumed constant 

IS pressure combustion process, if necessary, using the procedures established in 
PTC 4.1. 

In a similar manner, if the boiler absorption efficiency decreasestnon-stack 
losses iiicrease), consideration should be given to tenns affecting this cfTiciency : 
radiation & convection losses, gas turbine bearing losses if extracting power 

20 directly from the combustion gases, heat exchanger water/steam leaks, heat 
exchanger effectiveness, etc. The boiler absorption ^efficiency also has the 
potential to be determined with high accuracy. Generally this term approaches 
unity; thus its error is no greater than its compliment (if t}^ = 98 percent, one may 
safely assume a maximum error of ±2 percent). Although is dependent 

25 (through the term ERC) on Tie; a given degradation in Tic effect Ti^, the 
impact on relative changes is generally small. Also, by an iteration technique, ti^^ 
may be resolved without a priori knowledge of fijel flow rate. Thus, both t^ aiKl 
t|a, therefore T|boiieo i^ay be determined independent of fiiel flow. 

The enthalpy of the products (HPR) may be accurately calculaied usiiig 

30 methods discussed in the ^470 patent. The enthalpy of the reactants (HRX) may 
be accurately calculated using methods discussed in the '470 patent, provided the 
energy associated with injected limestone is included, if applicable, using similar 
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techniques. 

The Preferred Embodiment employs water and steam properties which are 
internally consistent within ±0.0008% maximum for enthalpy t:alculations. They 
derive from the 1969 work by Keenan, Keyes, Hill and Moore, modified by 
S NASA. Combustion and atmospheric ^ases have total consistency <zero error) 
given their explicit solution techniques between heat capacity, enthalpy and 
entropy; they are based on American Petroleum Institute standards. Moist air 
psychometrics are consistent within ±0.0010% for common enthalpy calculations. 
Below 38F, as low humidities are approached, consistency is within ±0.0090%. 

10 Heats of Formation are obtained from the world-recognized CODATA standards 
for thermodynamics. 

The above description of system equations and related procedures results 
in the concept that fuel flow rate may be determined with complete consistency 
from thermal efRciency. The accuracy of the computed fuel flow is intrinsically 

IS dependent on fuel chemistry and heating value, addressed below. The 
computation of boiler efiQciency using any operational and consistent method 
available is acceptable, provided it is independent of fiiel flow such that fuel flow 
may be back-calculated from the classical definition of thermal efficiency, and 
thus is thermddynamically consistent with Useful Enei;gy Flow Delivered from the 

20 system (and/or shaft power delivered from a combustion turbine). As used by the 
Input/Loss Method, the determination of boiler efficiency requtpes only 
computational consistency, not absolute accuracy. Computationally, boiler 
efficiency need only produce an accurate fuel flow, consistent with energy flows 
delivered BBTC (Eq.(5)). 

25 Combustion Equation and Fuel Chemistry 

Whatever specific procedures are employed to determine boiler efficiency 
(or system efficiency for a combustion turbine), this invention is based on a 
consistent treatment of the combustion equation. Any determination of efficiency 
is adequate vAien applied for this invention, provided the terms of the combustion 

30 equation are treated in an operational and consistent manner. The following 
combustion equation is taken from the EX-FOSS computer program for 
illustration of the preferred procedures. However, another consistent and 
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available method for the determination of efficiency may be found in Chapter 9 
of Steam, Its Generation and Use . 40th Edition, edited by S.C. Stultz and J.B. 
Kitto, published by the Babcock & Wilcox Company^ Barberton, OH. As applied 
in the present invention, ''consistency" is defined by three attributes: 1) use of 
S accurate molecular weights and consistently used thermodynamic properties; 2) 
use of accurate and correct molar balances, proof of which is demonstrated when 
such molar balances are used as a bases of system mass flow balances yielding 
minimum error; and 3) the incorporation of these same molar quantities within the 
computation of boiler efliciency. This later attribute is most important to the 

10 present invention. The concern is not the denominator of the boiler efficiency 
equation and whether boiler credits are included or not, but rather the ^equation's 
numerator. The terms HPR and HRX discussed above, as integrally based on the 
molar quantities defined below, with HNSL, defmed the uiseful power and/or 
energy flow delivered by the system per unit of fuel; thus use of these terms 

15 demonstrate required consistency. The equality mAp(HPR - HRX - HNSL) = 
BBTC is part of this consistency. 

In the present approach, the As-Fired fuel chemistry is first calculated 
based upon system measurements and the use of a mathematical description or 
model of the thermal system. All calculated system parameters (e.g., thermal 

20 efliciency, total effluent flow, emission rates, fuel flow, heating values and 
parameters associated with Second Law analysis) are thereafter calculated 
principally from the As-Fired fuel chemistry, and are therefore necessarily 
"consistent" with each other. By contrast, in each prior approach the various 
system parameters are either not calculated at all, or -calculated using a variety of 

25 calculational techniques and with different input information, with the result that 
the system parameters typically are inconsistent with each other to^ome degree. 

The preferred stoichiometric equation relatii^ reactants to products is 
presented as Eq.(29). The nomenclature used in Eq.(29) is imique in that brackets 
are used for clarity: for example, the expression "a2[H20]'* means the fuel moles 

30 of water, algebraically simply a^i the expression "d[C02]" means the -effluent 
moles of C02» algebraically simply ''d*\ Note that molecular quantities are used 
versus atomic for nitrogen, oxygen, and hydrogen; conversion is routine. The 
quantities comprising the combustion equation are traditionally based on an 



wo 99/49399 



PCrAJS99/a4S65 



-34- 

assumed 1 00 moles of dry gaseous product. This assumption is useful when 
measuring stack emissions since the commonly measured concentrations are based 
on dry molar firactions; if measured on a wet base, such measurements may be 
converted to a dry base. 

5 x{ao[CvRH2Rj + a, [NJ + ^^[HjO] + a3[0 J + a,{C] + ttsfH^J + a,[S] 

+ a,[COJ + a8[CO] + a,[HjS] + a,o[ash]U.pi^ 
+ MH201,„4^+ {(1 + P)(a[OJ+ a*A«[NJ + b^PH^ODU 

PLS 

= d[COJ + g[OJ + h[NJ + jP^O] + k(SOJ 
10 + {e[CO] + flHJ + 1[S03} + m[NO] + pLNjO] 

+ q[NOJ + t[Cyp,H2p,] + UCCYwHzpjlMinorComponcnB 

+ abpLs[CaS04-2H20] + {(1 - a + Y)bpLslCaO]}E,^pLs 
+ xa,o[ash] + w[CRei«J v{CRefuJ 

+ {pa(0J + Pa<|>^d + Pb^[H20lU,^ (29) 

1 S The following is presented to assist in understanding Eq.(29)» where all are molar 
quantities unless otherwise indicated: 

X — As-Fired fuel per 100 moles of dry-^as product (the solution base), 
tti = Fuel constituents per mole of As-Fired luel (wet-base), Eaj =1.0. 

ttj are derived from weight fractions as input to FUEL.EXE. 
20 b^ Moisture in entering combustion air, per base. 

bx s (o(1.0 + 4>Act)aN;vi/N„20f where a> is specific humidity of the local 
atmosphere determined from psychrometric measurements. 
Pb^ = Moisture entering with air leakage, f>er base, 
bz = Water/steam iri-leakage from the working fluid, per base. 

25 bz =xNAFE[m|eakage/(mAFNH2o)] 

bpLs = Molar fraction of pure limestone [CaCOj] 

required for zero CaO production, per base. 
bpLs = kp - k^c,; where kp s x(a6-a9)/ [L0+<l/kF)i„p„J; where the 

ratio of SO3 to fuel generated'S02, i.e., before reaction with 
30 limestone, (l/kp)j^p is supplied input, is the measured 

effluent ^Oj; without limestone injected, kp - k^ci. 
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Y " Molar ratio of excess CaCO, to stoichiometric 

CaCOj (e.g., 7 = 0.0 if no effluent CaO). 

Y = (n^LS^mAF)x N;^F / bpLs Ncacoa) - 1 0; where m^s is the system's 

indicated actual limestone mass flow rate, and ^ is a 
5 mass ratio of actual limestone to pure CaCOs it contains, 

z = Moles of H2O per CaS04, supplied as input based on 

periodiclaboratoiy analysis of boiler refuse. 
0 - Kronecker function: unity if (o^ + a,) > 0.0, zero if no sulfur 
is present in the fuel. 
10 jJ = Air pre-heater dilution factor a ratio of air leakage to tn» 
combustion air, a molar ratio 
P M/?act-10)/[a/?Act(1.0 + (|)^et)]. 

R^a = Ratio of total moles of dry gas from the -combustion process before 
entering the air heater, to the total moles of dry effluent (after 
15 leaving the air pre-heater); the air pre-heater leakage factor 

typically obtained fromieither the ratio of COj or Oj across the 
air pre-heater. 
a = Molar fraction of combustion Oj, per base. 
<j)Aa = Ratio of non-oxygen gases (nitrogen and argon) to oxygen in the 
20 combustion air local to the system. 

^Bci = Concentration of O2 in the combustion air local to the system, 
n^ = Molar quantities of drv^as product at the system's boundary for 
specific compounds: d, e, f, g, h, k, 1, m, p, q, t and u. 
25 The sum is denoted as En^. Dry airgases associated with 

system leakage (pa and pacl)^^) arc treated separately. 
= Molar quantities of non-gas product at the system's boundary for 
specific compounds: j, obpLs, (1 - a + Y)bpLs, xa,o, w and v. 
The sum denoted as En^. Moisture firom air leakage (pb^) 
30 is treated separately. 

Nj = Molecular weight of compound or element j, defined by the 
1 983 lUPAC to six significant figures. 
= Molecular weight of the As-Fired (wet) fuel. 
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input 



Ratio of mass flows of combustion air to As-fired fuel, based on 



V 



w 



indicated, uncalibrated, plant data, unitless. 
[(^''^Bottom.AJCCarbon/Bottom-Ash) 



5 



+ (f^^ny^(Carbon/Fly-Ash)] xN^f/Nc; WFi are weight fractions. 



Resolution of Eq.(29) proceeds by solving for all nj and nj, quantities. 
Minor component terms of Eq.(29) are typically resolved either through direct 
measurement (e.g.» for CO), or assumed zero values, or through assumed 
relationships (e.g., SO, is a function of an input molar ratio of SO, to SOi. defined 

10 above. All Minor Components typically have only low parts-per-million 
concentrations thus have little impact, with the possible exception of the terms v 
& w. The true importance and functionality of Eq.<29) to the Input/Loss Method 
lies in the fact that consistency of molar balances is needed Tor successful system 
understanding as discussed above, for conservation of mass flows, and for 

1 S resolution of fuel chemistry. For clarity the follovong major terms are associated 
with the system's thermodynamics: 



Eq.(29) is unique in describing at least three features of critical importance 
when improving the performance of conventional power plants with the Method 
of this invention. These include: 1) its abiliQ^ to address air pre-heater leakage 
30 through application of through the dilution term P; 2) the ability to desoibe 
effluent concentrations on either side of the air pre-heater, again through 



25 



20 



Total effluent water s J =j + b^p 

Total ef&uent oxygen s G = g + ap 

Total effluent = + En^i + P(a + a4>;^« + b^) 

Total effluents before the air pre-heater -Rj^n, -^R^JLn,,. 

Total As-Fired pure limestone = (1 + Y)bpLs 

Dry combustion air without air pre-heater leakage = (a + sl^^^^^ 

Dry air from air pre-heater leak^e present in effluent = P(a + a<t>Act) 

Wet combustion air = (a + ai^^ct + byj 

Total wet combustion air and air pre-heater leakage 



= (l + P)(a + a<l)^ + bJ, 
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application ofRAei 3) its ability to address limestone injection commonly used in 
fluidized bed combnstors; and 4) the use of an explicit term allowing for 
variable Oj concentration in the system*s local combustion air. Air pre-heater 
leakage dilutes all combustion effluents with moist air from the local environment* 
5 thus all important effluents HjO, CO^ and used for this invention are altered 
Furthermore, many times a power plant's more precise efSuent measurements, 
especially O^, may be found on the air pre-heate-'s inlet (economizer outlet), and 
not at the air heater outlet; refer to FIG* 6. Although most environmental 
regulation requires effluent measurements at the system's boundary^ translation 

10 between the air heater inlet and outlet measurements is many times essential. 
Eq.(29) allows for such translation through the term, defined above sudi that 
1 00 moles of dry gas are computed both at the upstream and downstream locations 
of the air pre-heater; see "Boiler** of FIG. 6. Thus effluents may be used by the 
present invention either upstream or downstream of the air pre-heater, refer to the 

15 G and J terms defined above, allovvdng conversion between measurements with 
and without air leakage. Combustion gas conditions upstream of the air pre-heater 
and before the exit of the heat exchangers and combustion region, see FIG. 6, 
would employ terms of g/J^ct and ji?Act Combustiongascs downsUream of the air 
pre-heater typically exit the system to the environment, in other words thegaseous 

20 effluentboundaiy of the system. Limestone injected into the combustion process 
will create additional effluent COj, and could decrease the-effluent H2O if the 
sul&te product is matrixed with water, CaS04 zH20. Of course CO2 and H2O «e 
two important effluents to the present invention. In addition to the basic 
stoichiometrics associated with limestone injection, the detail afforded by Eq.(29) 

25 allows numerous and obvious determinations of molar and mass ratios. For 
example the molar ratio of fuel sulfur to limestone*s<calcium, termed the limestone 
conveision ratio, important to a fluidized bed operator, is simply VK i Y)t>PLsl* 
Ihc concentration of Oj in the system's combustion air (termed A^, leading to 
the term affects all system effluent and combustion air terms. 

30 In summary, the aforementioned technique describes the process of 

calculating boiler efficiency of both any thermal system and in particular of 
practical conventional power plants based on effluent measurement data, fuel 
heating value and several parameters of minor importance. The next stage of the 
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process involves the recognition that a given fiiel has an iinique chemical 
composition, thus when burned will yield unique concentralions in its gaseous 
effluent. The principal gaseous effluents from any fossil combustion process are 
Nj, CO2, H2O, O2 and SOj. H^O, when effluent from combustion, is in its super- 
5 heated phase, thus acting as a gas. The source of is principally the air used to 
bum the fiiel and it has little chemical reactiveneiss, thus hs sensitivity to the iuel's 
chemical composition is not significant SO2 effluent concentrations are generally 
in the parts per million, thus its impact has minor importance. 

The relative concentrations of carbon, a^, and hydrogen, a^, foimd in any 

1 0 fossil fuel will have significant impact on the relative concentrations of COj and 
HjO found in the effluent. In addition, these effluents will be influenced by the 
following: O2 used to bum the fuel (i.e., the Air/Fuel ratio); fuel water, in- 
leakage of water; and water in the combustion air. This implies that the molar 
fractions of COj, HjO and Oj present in the effluent (the system's boundary, i.e., 

IS its smoke stack or translated from data from the air pre-heater inlet) must be 
unique relative to the fiiel input and supplied combustion air. 

The following elementaiy molar balances may be derived from the 
combustion Eq.(29). The Ff^ expressions are simply convenient groupings of 
quantities, principally comprising measured effluents (known values) which have 

20 the greatest influence on the individual -element of interest. Minor fuel terms 
carried within expressions are multiplied, initially, by an estimated fuel moles, 
X. These minor terms are quickly resolved when converging on x. An exception 
is r^, discussed below. Given these groupings the expressions of Eqs.(36) 
thru (40), with solution of the term "a" as discussed below, may be treated as 

25 known quantities. The elementary wet fuel components typically associated with 
coal are considered unknowns, as are the fiiel moles, these include the following: 
oi^y tts, oc^, a,o and x. Many coal-lired and bio-mass-fired plants use 
supplementary firing with gaseous fuel such as methane; others mix coals, 
frequently to adjust composite sulfur content; others mix hydrocarbon fuels such 

30 as coal, shredded tires, trash, and the like. The Input/Loss Method may account 
for all non*coal type hydrocarbon fuels, but they are assumed to have known 
chemistries with known fuel flows. Mixed fuels (such as multiple coals) are input 
to the system simulator (e.g., EX-FOSS.EXE) as a composite, combining them by 
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using the FIJEL.EXE program, although the resultant computed fuel chemistry, 
after proper weighting, is assigned to one fuel (e.g., the most variable coal), by 
assuming the other fuels have fixed or predictable, and known, chemistries. 
Alternatively, such assignment may be made to a combined subset of fuels. 



5 xa, =rNi-a4>Ac. (30) 

x(aj + as) =r„o (31) 

x(o3 + Oj/2) =ro, t32) 

xc^. =rcoj (33) 

xcte »rso2 (34) 

10 xa,o *-r;rt-x(«i + ai + 2aj + as/2- 1.0) (35) 

where: = 1-0 - (dAc,+ 6^^+ f + Gyvc,+ k^c, + 1 + tn/2 + q/2 + t + u) 

-<|)Ac.(/?Aa-10)/[/?^cl(4>Ac.+ 1.0)] (36) 

r,oo =(J^- b^P) + f+ZPAt/2 + ZP;,u/2.bz-bA 

- x(ZRO(»^ + a,) + bpLsOZ (37) 

15 r« = d^ + e^ + iG^-m + iJ«.-Km 

+ kA«+ 31/2 + m/2 + p/2 + q + 2r - a - b^^ - ba/2 

+ (3o - 2 r 2y + oz)bpLs/2 - x(a7 + ag^) (38) 

^coz = dA« + Cac + YP^t + YPbU + V + w - (1.0 + Y)bpi5 

- x(YRao + a, + ttg) (39) 

20 Tsoj = k^c, + 1 + YbpLs " xa» (40) 



=a-(GA«-ap) + Cac,/^ + fi^ + w + V - bpi5(^ - 1/2 - m/2 
- p/2 - q - 2r + x[Oo(YR - ZIWl) + a, + a,Q + Oj/Q] 

<41) 

In these relationships the subscript "Act" means an effluent measwsement or 
25 assumption (an "actual" value). The term in Eqs.(37) and <38) relating to the 
moles of effluent HjO (or jAe/^Act refonng to air pre-heater inlet cmditions). 
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could be input as a constant value or measured, ^iven the option su^ested by 
error analysis. All other values are either evaluated explicitly based on input data, 
internal models and/or have minor import but carried in the formulations for 
consistency. Eq.(35) is numerically correct, however if it is reduced by 
5 substituting x with Xmaf, and with a^AF^ see Eqs.(50) and (5 1 ), it is seen that 
Eq.(3 5) offers a trivial solution for aM^F-ich this result illustrates the null effects ash 
has as a neutral substance, and thus its resolution difficulties. 

As a group, these relationships are of critical importance for understandix^ 
the Input/Loss Method. If fiiel chemistry is to be resolved, thus heating value and 

10 thus accurate fuel flow, etc., then stoichiometric relationships ^generally 
representing Eqs.(30) to (41) must be resolved. These equations ar« not unique 
in their grouping of terms, further reductions and/or complexities are possible. 
The grouping of terms adopted here is principally the right-side terms of Eq.(29). 
For the following discussion, assume that fuel ash has zero^concentration. 

15 With this assumption, Eqs. (30) through (34) yield five equations with eight 
unknowns. Unknowns include a, through oc^ the terms "a" and "x'\ The term "x" 
is a convenience term and could be divided through changing the base of Eq.(29) 
to unity moles of fuel, thus eliminating use of terms comprising two unknowns, 
X and any aj. The requirement J^j^^jAj = 1.00, assuming ash is temporally 

20 ignored, is again a convenience; again the entire Eq.(29) could be divided through 
by or xa4, setting a base of unit moles of fuel carbon. Effluent N2-could be 
resolved by difference assuming 100 moles of gaseous effluent XCO^, HjO, Oj, 
SOj, and minor pollutants are measured or assumed), or Nj could be measured 
directly. By assuming a constant value for fuel nitrogen, a„ a minor fuel 

25 constituent, the unknowns are reduced to five with five equations. However, close 
examination will quickly indicated that this mathematical system invokes a stiff 
matrix. If ttj, again a minor fuel constituent, is assumed constant then Eqs.(31) 
and (32) represent two equations with two imknowns, and Uy Eq:(33) could 
be solve outright resolving its one unknown These manipulations are 

30 discussed to emphasize that algebrsdc manipulations must address the physical 
reality of the thermal system as done by the Input/Loss Method. The resolution 
of fuel carbon and fiiel hydrogen, of critical importance in resolving tuel 
chemistry, should address their intrinsic chemical relationship as a hydrocarbon, 
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and the influence played by the principal effluent measurements, d^^p and G^et- 
Invoking such relationships reduces the Method's sensitivity to instrumentation 
errors. The Preferred Embodiment employs this approach, relying on the belief 
that coals within a generic category (Rank) and geographical region, have well 
5 defined MAF carbon to MAF hydrogen relationships. Further, the fuel mole term, 
X, is used in the Preferred Embodiment for this single variable appears in all 
stoichiometric conversions to mass flows (the consistent determination and 
conservation of which is an object of this invention); and represents a iterative 
parameter for Eq.(29). Further, the Preferred Embodiment does not require that 

10 the minor fuel constituents be assumed constant. Further, as will become 
apparent, the Preferred Embodiment allows use of multidimensional minimization 
techniques which addresses instrumentation errors. 

Returning to the problem as posed by Eq.(29), including fuel ash, if fiiel 
chemistry is to be resolved then fundamental problems require solution - 

1 5 independent of algebraic manipulations. These problems include the following: 
1) for every unknown fuel quantity, a^, there is always, at least, additional 
unknowns in the combustion air term and the fuel moles, "a" and "x"; and 2) the 
fuel's solid non-combustible minerals (commonly referred to as "ash*') may not 
be resolved through effluent measurements since it is a neutral non-gaseous 

20 substance, the same input of xa,o in the &el appears in the effluent, further, its 
outlet flow from the system is not possible to measure directly with any reasonable 
accuracy. 

The first of these problems is solved by reducing Oj quantities to a molar 
MAF bases, eliminating the influence of the two components not chemically 

25 involved in the hydrocarbon fuel - water and ash - and then solving for a^^p.^ 
terms algebraically employing correlations of molar MAF fhictions of fuel 
hydrogen, nitrogen and oxygen as a function of carbon. Fuel MAF carbon is 
resolved explicitly by first solving Eq.(33) for x, converting to a MAF base, 
subtracting Eq.(32) from Eq.(3 1) to eliminate the terms, and into the 

30 result substituting Eq.(33)'s solution for Xmaf> yielding: 

ttMAF^ =-rco2(A5-2A3)/[rco2(B5-2B3)-r„20 + 2ro2] (42) 
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Note that the w & v terms of Eq.(39), related to reject & refuse carbon may be 
expanded, as functions of a^^p^ (see Eq.(29) notes), if their magnitudes aflBsct the 
accuracy of aj^p^. Given a^AF^, is resolved via Eq. (56), based on Eq. (33). 
In Eq.(42) the constants Aj and Bj are developed from reference fuel data 
5 associated with coal samples taken from the actual As-Fired, or obtained from 
generic specification of the coal Rank: 

^MAF^l ^A, +B,aMAF-4 <43) 
C^MAFO = Aj + BattMAF^ (44) 
ttMAF-S = As + BsttMAF^ (45) 

10 Fuel MAF sulfur is resolved explicitly by solying Eq:(34): 

^MAF-6 ~ ^502 / XjviAF (46) 

Fuel MAF water is resolved explicitly by solving Eq.(3 1) for x, and substituting 
into Eq.(32), yielding: 

.aMAF-2 = (aMAF-aTreo - aMAF-sTcj) / (Toj - (47) 

15 Both fuel MAF nitrogen and oxygen concentrations are small, and 

typically may be fixed as constants, or described by Eqs.(43) and (44). The MAF 
sulfur concentration may also be computed assuming 0^^^^^ = A^ + B^a^^ji, if the 
effluent SO^ measurement is questionable (or if limestone is injected and its 
conversion rate is questionable). As a further alternative, -either the oxygen or 

20 sulfur could be used to assure that Ej«o.i,3..9aMAF-j = I 0 (assuming Cmaf^ ^maf t^ 
<^MAF^ <Wf^ c>r any other secondary fuel has zero or known concentration). The 
sole criteria in deciding the exact methodology is the reliability and availability 
of effluent data and its relative impact on MAF fuel terms. For example, nitrogen 
is a major effluent, however using Eq.(30) to solve for fuel nitrogen, typically a 

25 very minor component of the fuel's makaip, would invite even slight errors made 
in determining effluent (made either by direct measurement or determined by 
difference suggested in Eq.(36)) to greatly amplify the uncertainty in Om^f.,. 
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Indeed, the F^j ^^"^ used to resolve the combustion O2, the term "a" via Eq. 
(58); wherein the term Xh^fOmaf-i is resolved through Eq.(43)such that ro2, which 
is dependent on ''a'\ may be resolved. Besides carbon, an important term is the 
fuel MAP hydrogen concentration which must represent a valid functionality with 
5 carbon. Experience has shown valid fimctionality may be easily achieved since 
the hydrogenxarbon chemical bond is predominate with carbonrcarbon; indeed 
the variability of a specific coal Rank lies not in its base chemistry but rather in 
its water and ash. FIG. 3 is a sample plot of MAP hydrogen versus MAF carbon 
for a coal fuel having highly variable As-Fired heating value, water and ash 

10 contents. In FIG. 3, both axes are MAF molar fractions per mole of MAF fuel. 

To this point in the process, resolution of fuel chemistry has been possible 
by assuming a Moisture-Ash-Free (MAF) basis. In carrying this out, it is not 
critical to the invention that carbon be computed as the independent quantity. 
Hydrogen, or any other MAF fiiel element could be considered as the independent 

15 variable. For example, if solving for hydrogen first, then Umaf^ ^/Cc^maf s)- 
However, the preferred procedure is to place independence with carbon, where of 
course the greatest accuracy and sensitivity may be found in the effluent 
measurements. CO2 does not exist in the combustion air to any appreciable 
concentration, it does not leak into the system, it is generated only from 

20 combustion and as a major effluent has obvious sensitivity to resolution of tt^Ap^ 
ofEq.(42). 

The second problem, dealing with ash, requires examination of the total 
system. The only system effect of fuel ash is as a pure dilutive or concentrative 
influence on fuel, and of course on the heating value. From a qualitative 

25 viewpoint, as fuel ash increases at the expense of carbon (for example), the 
amount of combustion air required to produce the same effluent O2, given less 
carbon to combust, will decrease. Given a decreasing heating value (higher ash) 
increased fiiel flow is required to meet the same useful enei^gy flow. Thvs an ideal 
system parameter for such sensitivities, which is routinely measured at fossil- 

30 fueled systems, is the Air/Fuel ratio. Generally such sensitivities are reasonable, 
a 1 0 percent increase in ash for a common coal will cause a 5 percent decrease in 
the Air/Fuel ratio, although many system and fuel parameters may affect such 
sensitivities. The wet, mass base, combustion Air/Fuel ratio (termed AF)^ as 
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computed by the system simulator is developed as follows: 

AF = im^,, + mMoi«u J / m^F (48 A) 

+P)l(a + a4)^N^,, + b^N^„ol/(xN^F) <48B) 

Expanding the term xN^f in Eq.(48B): 

5 xNaf = x(EjH^9Njaj + N,oa,o) (49) 

and then employing the following definitions of MAF fuel moles and fiiel 
constituents: 

XM^FHx/(1.04-a^P,,+ aMAFoa) (50) 
aMAF.j ^ ai (l -0 + aMAp.2+aMAF-ia) (51) 

1 0 allows substitution of Eqs,(50) and (5 1 ) into Eq.£49) for x and Uj, cancelling the 
term (1.0 + aMAF-2 ^ma?^\q> as intended, and then substituting into Eq. (48B) 
yields a solvable form: 

(53) 

Simplifying Eq.(53) and solving for MAF fuel ash, c^af-io yields the following 
results. Note in Eq.(54) that a normalized Air/I^uel ratio is used, normalized to 
actual plant data, defined by Eq,(57). x^af is substituted using Eq.(36), 



«MAF.|0 - rAshaMAFV(rc02N,o) " ^j^^ .^flM^T./t^ 10 (54) 

20 where: ^ (10 + P)[(a + a<|>A„)NAi, + b^NH^o] /^^a« <55) 

^MAF ~ rco2 / ^MAF^ "(56) 

AF^a = AF^ {AF^nlAF^td (57) 

a = (Tnj . XMApaMAF-l) / 4>Act (58) 



xNaf - X^,AF(^J«0...9NjaMAF.j + NjoaMAP-lo) 

= (1 .0 -h P)((a + a<|)Act)NAir + bAN^ool 
/ [xMAF(^ro...9NjaMAF.j + N,oaMAF.io)] 
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The variable AF^^ is the wet Air/Fuel ratio from the system's data collection 
devices (an unsealed signal); AFf^f^ is the wet Air/Fuel ratio consistent with 
reference fuel chemistry; and AF^^^ is the wet Air/Fuel ratio consistent with 
reference fuel chemistry but as would be produced from the system's data 
5 collection devices (a signal used for scaling). The value of N,o in Eq. (54) is input 
as a constant, or fitted as a function of aMAF.io(^us solving a quadratic equation), 
or fitted as a function of HHY^ap- Note that a system's indicated fuel flow 
measurement could obviously be used in place of AF, applying similar techniques 
as demonstrated in determining ttMy^.,©. However, use of AF is preferred since it 

10 integrally involves effluent and combustion air terms (through Fcoz, T^i and T^, 
and thus through such dependencies allows error analysis techniques to be 
operational and practical. 

It is noteworthy that the explicit procedure of determining fuel ash, and 
through use of the term (1 .0 + a^p,2 + aMAF-io) of Eqs.(50) and <5 1 ), allows any 

15 errors made in fuel water, a^Ap j* to be off-set by fuel ash, otMAF-io- This must 
occur since any given quantity xCj (wet-base) must be equivalent to XMAfttMAF j 
(MAF-base); if not, such wet to MAF conversions would numerically cause 
inconsistencies in the computed Air/Fuel ratio. There may be only one AF - wet- 
base is the same as MAF. Further, the effect may easily be demonstrated through 

20 off-line study. 

The selection of specific Calculational options is dependent upon the 
characteristics of fossil fuels. As an example usirig coals of various Ranks and 
other fuels, FIG. 4 illustrates a table of options relating to fiiel water, fiiel ash and 
the fuel's typical variability. Note that these options are subject to application of 

25 the multidimensional minimization techniques described herein; for^xample an 
assumption of constant effluent H2O may be corrected using error analysis. 

Thus all fiiel constituents are therefore determined on a MAF bases. nGiven 
these values, the wet base fuel molar fractions are determined, the wet base moles 
of fuel (x), followed by wet base (As^Fired) weight fractions: 

30 ttj = ttMAF-j/ (1.0 + aMAF.2+ ^MA^.tO^ (59A) 

X =XMAF(10 + aMAF-2+aMAF.ia) <39B) 

WF^ = ttjNj / (EojNj) (60A) 
= fVF^ / (1.0 - fVF2) (60B) 
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The present Input/Loss Method has important advantages over the prior 
techniques discussed in the Background. The following table clarifies dififerences. 



Determinations of 
Fossil Fuel Chemistries 



5 




Water 




MAF Elements Comments 




Input-Output 


no 


no 


no 


v^nemisiry not meniionecL 




Heat-Loss 


no 


no 




v.,«ncmisuy assumeu. 




(Jutput/ix>ss 


no 


no 


no 


"Approximate" chemistry assumed. 




Buna 


no 


no 


no 


Chemistry assumed. 


10 


Fuel Water & 












Ash Meters 


yes 


yes 


no 


MAF chemistry assumed known, 
may adjust MAF heating value. 




Thermal 












Analysis 


yes 


yes 


yes 


Off-Line, laboratory use with high 


15 










accuracy 'gas measurements, 
with sgram samples. 




Munukutla 


yes 


ho 


yes 


No air leakage, fuel ash is assumed 
known. 




Pat 5367470 


yes 


no 


no 


Water iterated^ known dry 


20 










chemistry relationships. 




Pat 5790420 


yes 


no 


no 


Water iterated, known dry 
chemistry relationships. 




Input/Loss 


yes 


yes 


yes 


Explicitx:alciilation for water, ash 



and carbon, with error analysis. 



25 Heating Valyig Cqlculationg 

Variations to a reference heating value may be determined by calculating 
the differential effects of MAF reference versus actual MAF fuel chemistry. The 
As-Fired fuel chemistry could be used to produce the MAF or dry. The heating 
value calculation involves use of a correlation relating oxygen, carbon, hydrogen 

30 and sulfur contents to a MAF heatii^ value, correcting for ash to obtain a dry-base 
heating value, then finally correctir^ for water to obtain the As-Fired value. The 
As-Fired value is then corrected for a constant pressure process if appropriate. 
The preferred correlation is scientifically based on the chemical binding energies 
between hydrocarbon elements. The correlation is taken from the works of 

35 Ghamarian & Cambel, which is based in-part on the well known work of Szai^t 
and Szargut & Stryrylska. The references include: A. Ghamarian & A. B. 
Cambel, "Energy/Exergy Analysis of Pressurized Fluidized Bed Combustor**, 
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Proceedings of the Intersociety Energy Conversion Engineering Conference, 
August 8-12, 1982, pp. 323-327; A. Ghamarian & A. B. Cambel. "Exergy 
Analysis of Illinois No. 6 CoaP, Energy, Vol. 7, No. 6, 1982, pp. 483-488; J. 
Szargut, "International Progress in Second Law Analysis", Energy, Vol. 5, 1980, 
5 pp. 709-718; and J. Szargut & T. Stiyrylska, "Approximate Determination of the 
Exergy of Fuels", Brennstoff-Warme-kraft, Vol. 16, No. 12, December 1964, pp. 
589-596. The correlation is accurate to within ±0.7 p^cent .AHHVm^p^p 
deviation for over four dozen short- and long-chained hydrocarbon compounds. 
For hydrocarbons like coal, demonstrated below, normally having a low oxygen 

10 content, the correlation's accuracy is estimated at ±0.5 p^xsent. A similar 
correlation by the same authors exists for hydrocarbons with high oxygen content 
To further reduce these errors, the Input/Loss Method computes a term AHHV,^p. 

based on a reference MAF-base heating value, HHVMAF.Re<* and the reference 
concentrations of oxygen, carbon, hydrogen and sulfiir (established diiring 

1 5 initialization). The term HHVMAF-unwrr is computed via Eq.<52) based on actual 
MAF concentrations of Eqs.(42) thru (46), or otherwise determined as discussed. 
Using Eqs.(63), (64) and (65) the on-line As-Fired heating value, HHV^, is 
resultant. The term N^af molecular wei^t of the MAF-base fuel (without 
water and ash). 

20 AHHVMAF-<«clia=HHVMAF.Rcf 

- (- 178387.18aMAF.3+ 1 83591. 92aMAF^ + 78143.68aMAF.5 

+ 127692.00aMAFjRer/NMAF-Rcf (61) 

HHVmaf™ = (- 178387. 18aMAF.3+ 183591.92aMAF^+ 78143.68aMAF.5 

+ 127692.00aMAF JAOua,/ NMAF-Actual (62) 
25 HHVmaF = HHVMAF-uncon+AHHVMAFHJclta <63) 

HHV,^ « HHVMAr(1.0.Jf7V,o) (64) 
HHVaf = HHV^^l.O-JrFj) <65) 



The preferred correlations used to determine heating values for the pi^sent 
invention are based on chemical binding eneigies. "Studies have demonstrated that 
30 traditional correlations, such as the Mott-Spooner correlation based on Dulong's 
formula - well known in the industry - are not adequate. The Preferred 
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Embodiment of the present invention requires at least the coefficients used in 
determining heating value to fall within certain ranges associated with three 
principal constituents of coaL Studies have indicated that using the above 
preferred constants, which fall within the required ranges, reduces the standard 
5 deviation of five dozen wildly varying coal analyses from ±228 to ±92 ABtu/lbm. 
The ranges of these coefficients, i.e., multiples the molar fractions Oj in Eqs^61) 
and (62), or their equivalent wei^t fractions (for this presentation of raxiges, the 
symbol WT) represents percent weight of j), include the following: for carbon 
molar fraction 174160cc^yNft^ to 184970aortoi/N|^, or in weight percent 

1 0 carbon, l^SIVF^^ to 1 54PFF^^ for hydrogen molar fraction 7761 Ofi^,,ogcf/Nfi^ 
to 114910ahydiDgei/Nfud assuming molecular hydrogen, or in weight percent 
hydrogen, 385 JFFhy^,^ to STOPFFhydjogefi ^md for oxygen molar fraction - 
1 63 1 90aoxyge/Nfug, to - 1 82390aoxygpyNtud assuming molecular oxygen, or in weight 
percent oxygen, -51 WF^^^ to -51 WF^^^ These ranges are independent of the 

1 5 fiiel base, whether MAF, dry or As-Fired fuel constituents are used. 

Thus the As-Fired ftiel heating value has been detennined. If a correction 
for a constant pressure process (HHVP) is required to the computed HHV^p, this 
correction may be made using the procedure of PTC 4. 1 . Iterations of the above 
system, combustion,, thermal efficiency, fuel chemistry and heating value 

20 calculations are preformed repeatability converging fuel moles, x, the higher 
heating value, HHV;^,r, such that total fuel energy flow, As-Fired, is adequate to 
meet the Useful Energy Flow Delivered, and minor stoichiometric teims to within 
tight tolerances. 

Gross Error P^tertiop md Mass B^laxic^s 

25 After convergence, the Input/Loss Method employs the primary 

measurements dxcp ^Act AF^ (used to determine fiiel chemistry, thus heating 
value, etc.) in comparisons to the computed values, dcai, ^ca ™^ ^^cqi detennined 
from the system simulator, to deduce inconsistencies in measurements. However, 
it has been found that such comparisons indicate only gross inconsistencies, and 

30 not necessarily those associated with incorrect system parameters such as air pre- 
heater leakage (i^Act) or combustion air (^act)> ^cl when used in an on-line 
application were not of value. A more robust and operational error analysis 
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employs multidimensional minimization techniques, discussed below. 

It is an important objective of this invention to assure conservation of mass 
flows. Such conservation, using the temis developed by this invention, 
interconnecting boiler efficiency, shaft power and usefiiL energy flow and 
5 effluents, assures the theraial system operator that consistency of these 
computations is acceptable. Mass balances are specified in the following table, 
note that if inlet and outlet mass flows disagree by more than 0.2%, errors are 
considered significant. Such errors may be associated with computer 
programming mistakes and, most importantly, inconsistencies in the principle 
10 effluent measuremoits. 

Fuel Flow Rate, m;,P,Eq.(21) = BBTC/[Tii^,eXHHVP + HBC)] 
Combustion Dry Air Flow Rate = m^^l-^^Xa + a4);vct)NAi/(xNAF) 
Combustion Air Moisture Flow Rate= m^j:( 1 + P)bAN„2o/(xN^F) 
In-Leakage of Water and Steam = Tn^fbJ^i^Q/ixNj^p) 
1 5 Pure Limestone (CaCOj) Flow = m ^pf 1 y^b^^ ^^^^^^j/IxN^ ^ ) 

EINLET FLOWS 

Dry Gas Flow from the Combustion 
Region (i.e.. Boiler Effluent, 

assuming 100 mole base) = m^pl 0ONqJ{R^ xN^f) 
20 Dry Air Leakage Flow at Boundary = m^^a + a(|>Act)NAi/(xNAF) 

Combustion Moisture plus Air Leakage 

Moisture at Boundary = m^fO Pb^)NH2o^(xNAF) 
Reject & Refuse Carbon Flow = m^f(y + w)N^(xNaf) 
Calcium Sulfate with Water Flow = ^Af^^Fis^caso^.Tmc/i^AF^A?) 
25 Calcium Oxide Flow = m^K 1 • a + Y)bpLsNcao/(xNAF) 

Ash (Bottom, Fly & Dust) Flow ^Mk^i^^fiiAsiMA? 

EOUTLET FLOWS 



By knowing the fuel chemistiy and flow rate, and the complete 
30 stoichiometric relationships, calculating individual emission flow rates, nijp^jo-i 
(Ib/hr), may occur as follows: 
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where Oj is the molar fraction of an effluent species on a dry-basis. The term 9, 
derives directly from solutions or measurements of the right-hand terms of 
Eq.(29), for example *so2 = The term x is the moles of fiiel per mole of diy 
S gaseous effluent. The emission rate per species, in units of pounds per million Btu 
of fiiel energy input, termed ERj, is given by the following: 

ERi = lOHpecic^CniAFHHV^p) (70A) 
- 10*W/[xN^pHHV^pl (70B) 

The emissions rale may be evaluated independently of the As-Fired total fuel flow 
10 rate, Eq.(70B). However, the computational accuracy of the fuel flow rate, as 
determined using the present approach, intrinsically affects the emissions rate 
through HHVy^F, x and N^p. Further, the process described herein allows the 
determination of total dry volumetric flow, at standard conditions of^gaseous 
effluent, denoted by VF, as required by environmental regulations. VF is 
1 S determined by the following (in standard-fP/hr): 

VF = Pgjn^pNg^[xNAp] (71) 

where Pg^^ and Ng^s are the standard density and avers^e molecular weight of the 
effluent dry gas. Of course, to determine the mass flow of all-effluents Eq.(69) 
may be summed. 

20 Minimization of Instrumentation Errors 

If the aforementioned error checks exceed an error tolerance, or for routine 
analysis of instrumentation mors, the Method of ibis invention may employ 
multidimensional minimization techniques to adjust measured and/or assumed 
effluent data (d^a & ^Acd* measured and/or assumed system data (i.^., AF,^ 

25 R^^^ ^Aaf* resulting in revised fuel constituents. Although in the system of 
equations of the Preferred Embodimmt, the effluent is considered a known and 
an accurate measurement, having no discemable error, this is not required as will 
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be observed in the following discussion. The assumption of known effluent P2 is 
supported in practice since industrial Oj instrumentation, commonly employing 
zirconium oxide technology, is relatively inexpensive. Thus a given thermal 
system may install numerous O2 instruments, upstream and downstream of the air 
S pre-heater» to achieve a result free of discemabie error. Furthermore, the & 
parameters, important components of the combustion equation, may also be 
determined using O2 sensors. 

When using the techniques of the Preferred Embodiment for minimizing 
error in the air pre-heater leakage factor R^^, Oj measurements are required 

10 between the heat exchangers/combustion region and the air pre-heater, ^ee f ig.6, 
and preferably immediately downstream from any effluent ash removal equipment 
but before the air pre-heater. 

After converging on fuel moles and heating value, with a^iven effluent O2, 
the Method of this invention then utilities certain parameters which are unique to 

IS a given fuel to determine instrumentation errors. These parameters include the 
moisture-ash-free (MAF) heating value, and certain parameters (termed 
developed from Roughton's work of 1980. Study of laige coal -data bases has 
indicated a remarkable consistency in MAF heating value for a^iven coal of the 
same Rank and geogn^hical region; for example, a Pennsylvania sub-bituminous 

20 coal's HHVmaf P^r cent deviation was found to be less than 0.02%. Roughton 
found that the parameter was within ± 1 .0% with 95% confidence for "a wide 
range of power station coals". This has been generally confirmed, but again 
confined to individual coal Ranks and regions; variations of ±075% for unique 
coals have been found to be common. The parameters for the total Tuel, iuel 

25 water and fiiel ash are determined in the following manner: 

" JuwiNhjO " ^MAF-lheo/^MAF-K^ASh) ' (^Dry-theoiNDiy-FtiePHV q^) 

(72) 

Note that: ^V^f^O^vJ^Ory = ^V^aia^lAAF-FjJi^WiAF' 

30 Avaicr = Ju«,^^H2o/(XDry.theoJ>Io,y.Fue|HH^ (73) 

^Ash = XMAFW*MAF.I(iNAsl/(XD,y.lh«»ND,^^^ (74) 
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where Jihcor* dihe<m k.iheor^d are used to develop these terms as determined 
directly from the updated fuel chemistry, assuming theoretical combustion. Note 
that no air pre-heater leakage is assumed for theoretical combustion: 1 .0, 
p = 0.0, thus Jo«or = jcheor; and = 0.20948. 
S The preferred procedure then establishes correlations of the Li parameters 

as functions of the actual input effluent or system data. For example, the 
following correlation relates the computed fuel factor to effluent COj: 

L,^ = K„ -H K,A« + K,3di« (75) 

Eq.(75) is established by varying COj (d^J and calculating Lpu^, for three^cases, 
10 keeping all other effluent parameters constant. Such variation is done 
automatically so that the reference value, L^f.^^^ is encompassed within the 
correlation's range. K,|, K,2 & K,3 are associated fitting constants. In like 
manner, correlations are formed for other parameters: 

(^«er + jE,A5h) versus Jytei 

15 (^w^ct + ^aJ versus^^rt 

Lpuei versus air pre-heater leakage factor (the term /?Act) 
Lp^i versus O2 concenttation in the 

combustion air (v4aci)» where =1/(1 4>Act)- 

The Input/Loss Method determines all such correlations automatically including 
20 their proper ranges. Note that the parameter Lp^^^ was found to yield best 
sensitivities to efiEluent CO2, as opposed to Roughton-s evaluation of the same 
teim, Ipueb which is used for sensitivity with /?Act and Aj^, He defined Ip^,, based 
on dried fiiel, as the mass flow of dry gas produced from theoretical combustion, 
Xi>y.th««NDry<ias> divided by the total eneigy of the dry liiel, \y^.ti,J^i^fut^JHVory' 
25 For perfect effluent measurements, without error, this later form and that used in 
Eq.(72) are of course identical, they represent different sides of the same 
theoretical combustion equation using dry fuel Use of l^^ and Lf^^^ was found 
of importance to uniquely sensitize their respective terms. The Preferred 
Embodiment of this invention employs the summed parameters (Avaier ^ash) ^ 
30 a function of J^j^ and-^F^tp individually; versus Z^^^ = /(J^) and = f(AF^). 



wo 99/49399 PCTAJS99/0426S 

-53- 

This summation was found essentially constant for a given fuel, more so than its 
individual terms. The Preferred Embodiment of this invention also employs, in 
addition to the Lj parameters, the fuel's MAF heating value. Again, the MAF 
heating value has obvious import given its constant nature for a given coal, but use 
S of the dry and/or the As-Fired (wet) heating value could be applicable if sample 
coal was especially prepared for testing purposes leading to ^correction of 
instrumentation signals (corrections being applicable to normal operations using 
fuel with unknown chemistry and heating value). 

These terms then form individual error functions, F^, the minimization of 
10 their aggregate forms the Objective Function, solved using well knovm 
techniques. 

= [(L, - L^rd'L^f.A' • [ 1 .0 - Ko-i A3(HHV„ - HHVR.,J^/HHVLf J (76) 
Objective Function = Z) Fj (77) 

where L, are computed values based on actual data, L^^,, are corresponding 

IS reference values; and HHV„ are computed heating values, HHYk^^^ are 
corresponding reference values for MAF, dry and/or the As-Fired (i.e., n = 1,2 or 
3). With this formulation, the power plant operator has the option of choosing any 
one or all of the five effluent and system measurements, with up to three types of 
heating values each. Further, any other relevant effluent and/or system 

20 measurement could also be incorporated, such as effluent Oj [L^uc\ ~ / (gAct)]» 
effluent SOj [lso2 ^/Q^aJ]* and/or indicated fuel mass flow [L^j^ = / (my^p^g^,)]. 
It should be noted that this error analysis would not be possible without the EX- 
FOSS program's ability to compute consistent effluent data based on given inputs 
of fuel chemistry and system parametos. Measured effluents are not input to EX- 

25 FOSS, but only effect computed fiiel chemistry. Capitalizing on this, the Method 
of this invention then forces consistency with fuel chemistry by correcting 
effluents (and thus chemistry) and system parameters throi^ error analysis. The 
end result is consistent fiiel flow, effluent flows, emission rates, boiler efficiency 
and system efficiency. 

30 Note that the term "reference fuel characteristics" includes an average or 

typical fuel chemistry and associated MAF heating value, preferably based on 



wo 99/49399 



PCT/US99/04265 



-54- 

historical data collection of ultimate analyses. The term "elemeatary" when used 
herein, referes to a fuel's composition by individual elements ^denoted by the term 
a, in Eq.(29)). The values oTL^,^^ Acf-waicr and 4ef.Aih included as a portion 
of the reference fuel characteristics, computed using the reference fuel chemistry. 
S Reference fiiel characteristics also includes whether the variability of fuel water 
and fuel ash in the As-Fired condition is predictable, or not. For any^iven fuel: 
fuel water may be held constant (including zero); fuel ash may be held constant 
(including zero); an accurate functionality may be observed Tor eitfier or both 
[e-g-« 0Wf.io ~/(HHVmaf)1; and/or fuel water and/or fuel ash may be treated as 

10 unknowns, determined using procedures of this invention. All of these possible 
variations for the treatment of fuel water and ash are included as a portion of the 
Method. Reference fuel characteristics also contain reasonability limits of the 
computed elementary fuel constituents. Reference fuel characteristics also contain 
fitting constants (e.g., K,,, K,2 & Kjj) associated with all correlations iielatiiig 

1 S dependent quantities to system measurements. 

It was found that the Preferred Embodiment used for the minimization of 
the Objective Function is the Conjugate Gradient method, described in the book: 
W.H. Press, S.A. Teukolsky, W.T. Vettering & B.P. Flannerv. Numerical Recipes 
in FORTRAN 77. The Art of Scientific Computing, Cambridge University Press, 

20 Cambridge and New York ( 1 992). This technique is sufficiently robust as not to 
require side constraints such as limiting boundary conditions. It produces effluent 
and system parameters (e.g., COj. HjO, /?, etc.) which as a set of data represents 
minimum error, thus correction factors are determined: for the COj sigrial 
dcor/^Act^ for the indicated Air/Fuel ratio AF^^,^^ = AFq^, /AF^^^-ctc. This 

25 technique, as others, require the partial derivatives of Eq.(76), thus: dFJdA^ 
dFJdJ^ dF'/dAFj^ct> dFJdRj^ and dF/dAj^^^ are computed in a routine manner. 
Other published techniques are common, such as those described in the book: 
G.N. Vanderplaats. Numerical Optimization Techniques for Engineering Design. 
McGraw-Hill Book Company, New Yoric (1984). Given the Method's procedure 

30 of allowing for effluent and system measiirements to be compared to 
corresponding, but consistently computed values based on updated fiiel 
constituents, common minimization techniques involving constrained searches, 
could be alternatively applied. These include Sequential Linear Programming, 
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Direction Set using Powell's method. Simplex method. Downhill Simplex method, 
Simplex method with product form inverse, Quasi-Newton method, and others. 
Commercial products are also available, such as from Lindo Systems, Inc. of 
Chicago, III. 



lies with use of neural network technology (NN)* NN recognizes patterns in data, 
which if attributed to certain effects, makes corrections to the causal data. 
Specifically, data from instrumentation required for the Input/Loss Method xould 
be analyzed for effects on the fuel's computed results (e^g., functions), and 

10 corrections then applied such that thexomputed fuel chemistry yields consistent 
As-Fired heating values. Given such corrections to instrumentation signals, the 
established Method would proceed as herein described. Over checks couM be 
established which monitor a systems' fiiel energy flow and indicated fuel flow, 
comparing the computed with the measured. However, this is not preferable 

IS given the nature of NN technology as it is more amiable to large data bases, 
representing processes too complex for explicit thermodynamics. Numerous 
commercial NN software packages are available; for example from Neural Ware, 
Pittsburg, PA; California Scientific Software, Nevada City, CA; or The 
Math Works, Inc., Natick, MA. 

20 An example of using Conjugate Gradient minimization on a system of 

equations is presented in the following table using actual power plant data, 
minimizing the errors associated with d^^^ (and using HHVmap), Jasp ^acv ^aci 
(and using HHVd^) and Aj^ 



A further technique applicable to the reduction of instrumentation errors 



25 



Instrument Signal 
or Assmnption 
Dry stack COj 



Uncorrected 



J2ata 

12.619% 
6.300% 
6.980 
1.2500 



Corrected 



J2ata 

12.610% 
7.447% 
7.065 
1.2517 



Wet stack HjO 



30 



Assumed Air/Fuel 

APH Leakage (JJacO 
Assumed Oj/Air 



20.9480% 



20.8099% 



As an example of an alternative technique, <:onsider rkoui>^^i^"^^i^'^ 
affecting fuel carbon, water and ash as su^ested by Eqs.<37), (38), (39) and<55), 
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but based on the computed values dcab ^cai '^^ceh Obvioiisly the deviations 
between the measured and calculated values of a^,^p^ therefore may be 
determined: 

5 Any non-zero value of Eq.(78) may be ultimately attributed to inconsistencies in 
the efHuent and/or system measurements. Therefore a valid fuel composition is 
the one which minimizes these deviations thus being consistent with both 
stoichiometry and measurements. This is not preferred since no use is made of the 
Rj^ct '^Aci terms as might affect effluent measurements and/or independent 
10 parameters (e.g., Lpiiei & HHVmaf)- 

Svummm^ gf Thermodynamic Understandiine 

The present Input/Loss Method provides improved understanding of the 
thermodynamics of fossil systems, as compared with the prior approaches 
discussed in the Background. The present Input/Loss approach allows 
IS determination of complete As-Fired fuel chemistries, heating values, fuel flow, 
effluent flow, consistently calculated ash flow, -consistently -calculated emission 
rates, and consistently calculated system efficiency. Comparisons to other methods 
are presented in the following table: 

Complete As-Fired Calc. Calc. Calc. 

20 As-Fired Heating Fuel EfiHuent Ash Emission System 



25 



30 



Method Fuel Chem. 


Value 


Flow 


Flow 


Flow 


Rates 




Input-Output 


NA 


PM 


PM 


NP 


NP 


NP 


yes 


Heat-Loss 


PM 


PM 


NA 


NP 


NP 


NA 


yes 


Ouiput/Loss 


(no) 


PM 


yes 


PM 


NP 


yes 


yes 


Buna 


PM 


PM 


yes 


NA 


NA 


NA 


yes 


Fuel Water & 
















Ash Meters 


(no) 


(no) 


NA 


NA 


NA 


NA 


NA 


Thermal 
















Analysis 


Xno) 


(no) 


NA 


NA 


NA 


NA 


NA 


Munulcutla 


(no) 


(no) 


(no) 


NA 


NA 


NA 


yes 


Pat. 5367470 


(no) 


(no) 


(no) 


(no) 


NA 


yes 


yes 


Pat. 5790420 


(no) 


(no) 


(no) 


(no) 


NA 


yes 


yes 


Input/Loss 


yes 


yes 


yes 


yes 


y€s 


yes 


yes 
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Notes include: 

PM => Perfect measurement or perfect assumptions made. 
NA => Not Applicable, or not done given original pi^entation 
of method. 
5 MP =^> Not Possible. 

^^Calc." Calculated with consistency based on non-direct 
measurement or assumptions, 
(no) => Results arc possible but restricted to either known MAP fuel 
chemistry, and/or known fuel ash, and/or known fuel 
10 water; for Munukutla, his method is restricted toxonstant 

and perfect assumptions as to fuel ash. 



At the practical level of applying the Input/Loss Method at an operating 
power plant, its use of an air pre-heatcr leakage factor (Ra^df ^ pre-heater 
dilution factor (p), and a variable Oz/combustion air term (A^ thus <j)Aci)> which 

1 S allow use of effluent measurements on either ^ide of the air pre-heater, coupled 
with techniques whidi correct effluent instrumentation through terror analysis, 
assure operational and reliable themial performance monitoring. Considering that 
boiler efficiency is formed from fuel constituents (based on measured effluents), 
that a set of consistent effluents are computed based on these fueUconstituents, 

20 that the measured effluents are corrected to yield minimum differences with the 
computed through error analysis, the power plant operator has assurance of 
consistency. Thus this Method, considering its procedures, defii^ '^self 
consistency*'. 

When applying the Method of this invention for the improvement of 
25 thermal system efficiency the Preferred Embodiment is through an established 
technique using Fuel Consumption Indices (FCI). FCI are developed firomexergy 
and thermodynamic irreversibillQ^ computations of the "system^s major 
components and processes based on the Second Law of thermodynamics. These 
components and processes include in part: the combustion {Hx>Gess, combustion 
30 gas to wcricing fluid heat exchangers, the turbine cycle, and the process of direct 
generation of electricity. FCIs are computed by the EX-FOSS computer program. 
For the purposes of this invention this technique is defined as a **Second Law 
analysis of the thermal system'*. The most important inputs to this technique 
includes specification of As-Fired fuel chemistry, As-Fired fuel mass flow and 
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combustion gas mass flows; which given the consistent approach of the present 
invention, are fully resolved. Other required inputs include system mass flows 
including those of the working fluid, system enthalpy flow balances and 
miscellaneous data. Thus the Method of this invention are uniquely suited to 
5 allow the benefits of Second Law analysis of thermal systems to be implemented. 
The reference for this technique is: F.D. Lang and K.F. Horn, 'Tuel Consumption 
Index for Proper Monitoring of Power Plants", EPRI 1991 Heat Rate 
Improvement Conference, May 7-9, Scottsdale, AZ, also available from Exeigetic 
Systems, Inc. of San Rafael, CA in its latest version. 

10 Summary of Procedure 

In sununary, details of the procedure involves first several initializing 
procedures, followed by sequential procedures preformed by a computer operating 
in communication with the systems* data collection devices while it is operating, 
in real time. 

15 Prior to real time operation, the Input/Loss Method requires certain 

initializing data involving reference fuel chemistiy and heatii^ value, and 
stoichiometric data associated wth the reference ftiel chemistry. To supply such 
data a series of reference fuel chemistries, as ultimate analyses (identifying by 
weight fractions the individual elements), established from historical laboratory 

20 reports is desirable to determine relationships between an independently 
determined MAF molar carbon fraction aMAF-4> and dependent hydrogen, nitrogen, 
oxygen and sulfur fractions; i.e., ttMAF-j ^/C^maf J- ^ coal fuel, nitrogen, oxygen 
and sulfiir components are small relative to^bon, hydrogen, water and ash, and 
thus have minor importance to the fiiel's heating value. However, the relationship 

25 aMAF-5 "^/(^maf-a) should be established based on actual data. An associated 
reference higher heating value, consistently -established from the same historical 
laboratory reports acquired for the fitel, is desirable, reduced to a MAF base. 
However, if no such data are available, then the fliers minor components of 
nitrogen, oxygen and/or sulfiir may be held constant and/or calculated -explicitly 

30 if effluent data has the appropriate accuracy, while the relationship Omafo ~/ 
(ttMAF-4)» and the reference MAJF heating value may be established from coal 
analyses associated with the generic coal type. During the initialization, an -error 
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analysis is preferable using sampled data from the system to determine both the 
interaction of instrumentation errors on effluent COj. effluent HjO, Air/Fuel ratio, 
the term R^^ (and thus air pre-heater leakage), and the term and to confirm 
model selections related to fuel chemistry computations. 
S FIG. 1 illustrates the initialization processes. Box 1 10 indicates the initial 

(and one-time) setups for the indicated items. Box 120 depicts the process of 
developing a reference fiiel, involving all common fuels used by the system. 
Reference data is input as wet-base weight l&actions to a system simulator such as 
EX-FOSS.EXE, but also reduced to Moisture-Ash-Fiee (MAF) base for use by 

10 Eq.(61). Box 1 30 represents a simple organization of all initialization data. Box 
140 depicts using the same system simulator used when on-line, for off-line 
sensitivity studies resulting in selection of appropriate calculational options. 
FIG. 3 presents a typical relationship between MAF molar hydrogen and MAF 
molar carbon. This particular data is based on wildly varying coal fuel, whose ash 

1 S varied by sk 1 9 percent about the mean, whose fuel water varied by ±40 percent, 
and whose heat value typically varied by dbl2 percent All secondary fuels (e.g., 
coals and/or gas and/or oil fuels) are specified, and assumed 4o have constant 
values. FIG. 4 suggests Calculational options associated with various fuels, 
discussed below. Each of these options should be confirmed during the 

20 initialization process involving error analysis. Box 140 of FIG. 1. With these 
initialization steps completed, the analysis proceeds to real time, on-line, analysis. 

As used herein, the temi "obtained" or "obtaining** is defined as measuring, 
calculating, assimiing a non-zero value, assuming a zero value, estimating, 
gathering from a published reference work gathering from a data base, or any 

25 other operational approach 

Hiis paragraph discusses and defines '^operating parameters" obtained fi-om 
a thermal system applicable to this invention. . Effluent C02,<>2> NOx* and SO2 
measurements are made (at the system's effluent stream or **smoke stack"), or 
before the air heater (tenned the "boiler" side of the air pre-heater, i.e., upstream 

30 of the air pre-heater device). O2 measurements are preferred before the air pre- 
heater. Indicated wet-base combustion Air/Fuel mass ratio and effluent H2O 
measurements are required, or assumptions made, depending on the refeience tUel 
characteristics and appropriate error analyses. Measurements coniprising the 
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Air/Fuel ratio could be made on a volume base, or a dry-base, then converted. 
Effluent temperature is measured, that is the average temperature associated with 
the combustion gases at the boundary of the system (caution must be exercised in 
measuring non-stratified gas flows). The inlet/outlet ratio of CQzTpreferred), CO, 
S or O2 across the air pre-heater (these could be obtained off-line, based on periodic 
testing or judgement, but on-iine is prefeired), is used for the determination of air 
pre-heater leakage in accordance with ASME PTC 4.3, using the J^^ct^^^X^^ 
thus 0, see notes associated with Eq. (29)). Fuel temperature at an aj^ropriate 
thermodynamic boundary of the system is obtained. Air psychrometrics arc 

1 0 obtained at the boundary of the system (e^g., diy and wet bulb temperatures, or dry 
bulb and relative humidity). The dischaige temperatures of the air as it^xits^ach 
air heating or cooling device (but before it reacts with the fuel) are required; Tor 
example, such devices might include the air pre-heater, foiced-draA fan, steam-to- 
air heater, etc. Measurements are required to determine the total energy flow 

I S deposition to the working fluid fiom the combustion^gases. For a power plant, 
such measurements typically include Teedwater flow to the steam ^nerator, 
feedwater pressure and temperature, steam flow from the steam generator if 
different than the feedwater flow, steam pressure, steam temperature or quality (or 
assumed quality), and, if applicable, reheat flows, and reheat inlet and outlet 

20 pressures and temperatures. For a conventional power plant, determination of 
accurate reheat flows generally requires understanding of steam turbine flow 
distributions (involving shaft seals, feedwater heaters, bypass leakages, 
attemperation spray flows and the like). 

FIG. 2 presents the Calculational sequences associated with the Method 

25 when monitoring a fossil-fired system on-line, i.e., in essentially real time. Box 
250 represents data passed from the initialization process depicted in FIG. 1 . Box 
260 reduces fiiel data and prepares input for the system simulator. Reduction of 
fuel data involves combining the newly computed fuel with "fuels which have 
constant chemistries. Box 270 is system data as indicated; other minor data could 

30 be updated periodically while on-line. Box 280 is the system simulator which, 
given specification of a composite fuel, produces fuel flow and system-efficiency. 
Box 285 is the computation of fuel chemistry and fuel heating value for both the 
composite fuel, and given the fixed compositions of the secondary fuels, the 
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composition of the primary fuel. Box 287 tests for convergence of the process 
based on the composite fuel moles, heating value, and minor terms. In general, 
convergences are within 0.5x 1 0^ percent for the computed As-Fired fuel moles. 
Box 288 represents results from the error analysis indicating which instruments 
5 require servicing, and numerical minimization of instrumentation-errors associated 
with the input measurements of Box 270 versus those computed from Box 280. 
291 represents the transfer of correction factors produced from^error analysis for 
effluent instrumentation and/or data assumptions. Box 290 represents a decision 
as to whether the correction factors warrant update and^subsequent i^. Box -292 

10 determines effluent flow and emission rates, refer to EqsX^9) thru (71) and 
associated discussion. Box 294 provides the results of the Second Law analysis 
of the thermal system, other output and reports to system operators such that 
corrective actions may take place; also, reports are provided to regulatory 
authorities. Box 296 is a decision to turn the process off {quit) or not. Box 298 

IS is to quit. 

FIG. S presents a sample of calculational results associated with a power 
plant buming coal with wildly varying chemistry, as discussed earlier in relation 
to FIG. 3. This plant routinely uses coal which varies by ±12 percent in heating 
value. As may be seen in FIG. 5, the analysis produced average errors less than 

20 rtO.S percent in As-Fired heating value, a remarkably good result. 

FIG. 6 is a schematic representation of a conventional or fluidized bed 
power plant illustrating use of stoichiometric relationships important in applying 
this invention to practical power plants. It should be studied in conjunction with 
£q.(29) and its associated discussion. Specific discussion of FIG. 6 is presented 

25 earlier. Limestone injection is shown in FIG. 6 which is commonly used in 
fluidized bed power plants. 

In summary, by mass and energy balances based on unity fuel flow rate, 
by using accurate molecular weights, by using accurate thramodynamic properties 
of water, steam and combustion gases, by knowing the gross electrical and/or 

30 mechanical power produced from the combtastiongases directly or indirectly, by 
knowing the net energy flow supplied to the working fluid, and by recognizing the 
integral relationship of effluent CO2, HjO, SO2, O2, Air/Fuel ratio, air pre-heater 
leakage and O2. concentration in the local combustion air - to their fuel 
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constituents, revised fuel chemistry may be computed, the heating value based on 
the revised fuel chemistry may be computed, and fuel flow to the system may be 
computed. Knowing fuel flovir and heating value allows determination of affluent 
flows and emission rates. Knowing fuel flow, heating value and fiiel chemistry 
5 allows routine tracking of a fossil-fired plants' overall thermal efficiency, thus 
continuous correction of problems impacting thermal efficiency is possible. 

Although a particular embodiment of the invention has been described in 
detail for purposes of illustration, various modifications and<eiihancements may 
be made without departing from the spirit and scope of the invention. 
1 0 Accordingly, the invention is not to be limited except as by the appended -claims. 
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What is claimed is: 

1 . A method for quantifying the operation of a fossiI«fueled thermal 
system having a heat exchangers/combustion region producing combustion 
S products, the method comprising the steps of: 

before on-line operation, obtaining reference fiiel characteristics; and 
thereafter 

operating on-line, the step of operating on-line including the-steps of 

measuring a set of measurable operating parameters, including at 
1 0 least effluent concentrations of and COj, these measurements being made at a 
location downstream of the heat exchangers/combustion region of the thermal 
system, 

obtaining the concentrations of pollutants CO, SO2 and NOx, 
obtaining an effluent concentration of HjO, if reference fuel 
1 S characteristics indicate fuel water is not predictable, as an obtained -e'ffluent H2O, 
obtaining an indicated Air/Fuel ratio, if reference fuel 
characteristics indicate fuel ash is not predictable, as an obtained Air/Fuel ratio, 
obtaining the ambient concentration of O2, 
obtaining the air pre-heater leakage factor, and 
20 calculating the complete As-Fired fuel chonistry, including fuel 

water and fuel ash, as a function of the reference fuel characteristics, the set of 
measurable operating parameters, the'concentrations of pollutants CO, SO2 and 
NOx, ^he obtained effluent HjO, the obtained Air/Fuel ratio, the ambient 
concentration of O2, and the air pre-heater leakage factor. 

25 2. The method of claim I, wherein the st^p of operating on-line 

includes the additional step after calculating the complete As-Fired fuel chemistry, 
of 

calculating consistent moisture-ash-free, dry and As-Fired tuel heating 
values as a function of the complete As-Fired fuel chemistry and the reTerence 
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fuel characteristics. 

3. The method of claim 2, including, after the step of calculating 
consistent moisture-ash-free, dry and As-Fired heating values, the additional steps 
of 

S obtaining fuel factors associated with the As-Fired fuel chemistry, 

updating the reference fuel characteristics based on the calculated As-Fired 
fuel chemistry and its sensitivities to fuel characteristics, 

completing a multidimensional minimization analysis employing fiiel 
factors and heating values to minimize the collective error associated with at least 
10 one of the measured effluent CO2, the obtained effluent HjO, the indicated 
Air/Fuel ratio, the ambient concentraticxi of O2 and the air pr-e-heater leakage 
factor, 

calculating and applying for subsequent on-line analysis connection factors 
to the measured eflHuent CO2, obtained effluent HjO, the indicated Air^uel ratio, 
15 the ambient concentration of O2, and the air pre-heater leakage factor. 

4. The method of claim 2, wherein the set of measurable operating 
parameters includes effluent temperature, and wheFein the method includes an 
additional step, after the svsp of calculating consistent As-i^ired fuel heating value, 

• of 

20 calculating a consistent boiler efiBciency as a function of the^mplete As- 

Fired fuel chemistry, effluent temperature, effluent ^oncentratioris, and heating 
value. 

5. The method of claim 4, wherein the step of operating on-line 
includes an additional step of 

25 measuring the useful output derived directly from the combustion 

products, 

and the step of measuring a set of measurable operating pararneters includes the 
step of 

measuring a useful power and energy flow produced from a working fluid 
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heated by combustion products, 

determining an efficiency associated with the working fluid*s process, and 
further including an additional step, after the step of calculating a consistent boiler 
efficiency, of 

5 calculating a system thermal efficiency as a function of the efficiency 

associated with the working fluid's process and thexonsistent boiler efficiency. 

6. The method of claim 4, wherein the siep of operating on-line 
includes an additional step of 

measuring a useful output derived directly from the combustion products, 
1 0 and including an additional step, after the step of calculating complete As-Tired 
fuel chemistry, of 

calculating a consistent heating value as a function of the complete As- 
Fired fuel chemistry, and 

including an additional step, after the step of calculating a consistent beating 
15 value, of 

calculating a consistent As-Fired fuel flow rate as a function of the 
consistent tx>iler efficiency, the useful output, and the As-l^ired heating value. 

7. The method of claim 6, including an additional step, after the step 
of calculating a consistent As-Fired fuel flow rate and boiler^fficiency, of 

20 calculating a consistent total effluent flow rate as a function of the 

consistent As-Fired fuel flow rate. 

8. The method of claim 7, including an additional step, after the step 
of calculating a consistent total effluent flow rate, of 

calculating individual effluent flow rates as functions of the total effluent 
25 flow rate and effluent concentrations. 

9. The method of claim 8, including an additional step, after the step 
of calculating individual effluent flow rates, of 

calculating consistent individual emission rates as a function of the 
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individual effluent flow rates, consistent As-Fired fuel flow rate, and consistent 
As-Fired heating value. 

10. The method of claim 1, wherein the thermal system comprises a 
thermal system selected from the group consisting of a coal-burning power plants, 

5 oil-buming power plants, gas-fired power plants, biomass combustors, fluidizaed 
bed combustors, a conventional electric power plant, a steam generator, a package 
boiler, a combustion turbine, and a combustion turbine with a heat recovery 
boiler. 

1 1 . The method of claim 8, including an additional step, after the step 
10 of calculating consistent Individual effluent flow rates and those quantities on 

which they are consistency based, of 

calculating a Second Law analysis of the thermal system and Fuel 
Consumption Indices. 

12. The method of claim 1 1, including an additional step, aifter the step 
IS ofcalculating Fuel Consumption Indices, of 

adjusting operation of the system to improve its efficiency based upon the 

results. 

13. The method of claim 1, wherein the step of operating on-line 
includes the additional step aflor calculating the complete As-Fired fuel-chemistry, 

20 of 

obtaining an As-Fired fiiel chemistry including at least one fuel 
constituent selected from the group consisting of weight percent-carbon, weight 
percent hydrogen, and weight percent oxygen of the fliel; and 

calculating consistent moisture-ash-free, dry and As-Fired fuel 
25 heating values as a function of the complete As-Pired fiiel chemistry and the 
reference fuel characteristics in Btu per pound of fuel, the step of calculating 
including a step of forming products of numerical coefficients times the weight 
percent of the fuel constituent, wherein at least one of the numerical coefficients 
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is selected from the group consisting of from 145 to 154 for weight percent 
carbon, from 385 to 570 for weight percent hydrogen, and from -51 to -57 Tor 
weight percent oxygen, and combinations thereof 
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